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FOREWORD 


This  report  consists  of  seven  parts,  originally  written  so  that  each 
could  stand  alone  as  a  separate  article  or  technical  note.  Each  part 
contains  its  own  abstract,  which  is  placed  at  the  start  of  that  part. 
Likewise,  a  detailed  table  of  contents  is  placed  at  the  start,  and  ref¬ 
erences  at  the  end  of  each  part. 
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ABSTRACT 

X. 

From  measured  waveforms  it  has  been  determined  that 

the  electromagnetic  pulse  from  a  nuclear  detonation  is 

• 

associated  with  the  gradual  establishment  of  a  large 
dipole  moment,  which  remains  suspended  in  the  atmosphere 
at  the  conclusion  of  the  EM?  waveform.  From  the  gradual 
growth  and  the  large  magnitude,  it  is  inferred  that  the 
source-current  region  is  large,  and  that  retardation 
across  this  source -current  region  is  an  Important  aspect 
of  the  EMP  phenomenon.  Retardation  is  also  important 
in  Individual  processes,  since  the  primary  Compton 
electrons  are  relativistic,  and  the  .forward  directivity 
of  their  radiated  electromagnetic  fields  can  be  attributed 
to  retardation  and  can  be  encompassed  in  a  distrlbuted- 

, current  picture  which  allows  correctly  for  retardation,^. _ 

When  retardation  is  properly  incorporated,  it  is  found  that 
for  early  times  there  is  a  momentary  radial  electric  field 

-projected  above  the  detonation  center,  with  the  polarity . 

to  drive  secondary  electrons  upward,  thereby  contributing 
substantially  to  the  EMP  dipole  moment.  The  AFWL  non-caucal 
solution  gives  the  wrong  polarity  for  this  early-time 
radial  field,  but  the  non-causal  parts  of  the  AFV/L  program 
can  be  replaced  by  ©  causal  iteration  procedui-e ,  without 
affecting  the  parts  of  the  program  which  deal  with  electron 
attachment  and  air  chemistry,  and  with  the  gamma  rays  and 
the  primary  source  current. 
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1.  INTRODUCTION 

Previous  reports  In  this  series1*2*^  have  dealt 
with  three  aspects  of  the  EMP  (electromagnetic  pulse) 
problem.  There  is  first  the  phenomenological  aspect* 

From  measured  waveforms  It  can  be  determined  that  a 
nuclear  explosion  in  the  lower  atmosphere  establishes 
a  dipole  moment  in  the  air  with  a  negative  polar! by. 

The  detonation  raises  negative  charge  and  leaves  it 
suspended  in  the  air.  The  time  required  for  this  process 
to  take  place  is  on  the  order  of  100  microseconds, 
being  lonser  for  the  larger-yield  detonations,  shorter 
for  those  of  small  yield. 

The  growth  of  the  dipole  moment  is  gradual,  not 
abrupt.  Also,  detailed  analysis  of  the  measured  waveforms 
(Chapter  II  of  Referenoe  -2)  shows  that  the  negative  charge 
remains  suspended  at  the  end  of  the  EMP  waveform,  decaying 
to  zero  only  after  a  time  which  is  very  long  in  comparison 
with  the  duration  of  this  very-low-frequenoy  pulse.  These 
two  observations,  based  on  the  measured  waveforms,  lead  to 
the  inference  that  the  source  ourrents  are  not  mainly 
confined  to  the  near  vicinity  of  the  detonation  point,  but 
are  distributed  over  distances  measured  in  kilometers, 
at  least  in  the  vertical  direction  above  the  detonation 


oenter. 


The  seoond  aspect  of  the  EMP  problem  which  has  been 
dealt  with  In  detail  is  the  relativistic  nature  of  the 
primary  Compton  currents  which  provide  the  driving 
source  for  the  main  EMP  phenomena.  Gamma  rays  from  the 
nuclear  detonation  produce  Compton  electrons  through 
collisions  with  air  molecules.  Because  of  the  blockage 
by  the  ground  under  the  detonation  point,  most  of  these 
Compton  electrons  are  directed  upward,  and  they  thus 
contribute  to  the  establishment  of  a  negative  dipole 
moment.  At  the  same  time,  their  motion  (including  their 
initial  acceleration  at  the  moment  of  ejection  from  an 
air  molecule,  and  their  more  gradual  deceleration  &'s-  they 
are  slowed  by  ionizing  collisions  with  air  moleoules) 
will  produce  electric  and  magnetic  fields  which  can 
aot  on  other  electrons  in  the  violnlty.  These  other 
electrons  lnolude  not  only  the  primary  Compton  electrons 

*  i 

but  aiso  secondary  electrons  produced  by  the  ionizing 

collisions  which  slow  the  primary  Compton  electrons, 

and  other  secondary  electrons  released  through  the 

X-ray , 

photoionizing  action  of/ visible,  and  ultraviolet  radiation 
emitted  by  the  nuclear  detonation. 

Because  of  the  relativistic  velocities  of  the 
Compton  electrons,  there  is  a  forward  directivity  in  the 
electromagnetic  fields  that  they  generate.  This  forward 
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directivity  is  incorporated  in  the  field  expressions 

obtained  from  the  Li^nard-Wieohert  potentials/* 

However,  these  potentials  were  in  turn  derived  from 

• 

a  distributed-charge  pioture,  in  which  the  source 
electron  was  assumed  to  be  a  compact  distribution  of 
moving  eleotrlcal  charge.  The  retardation  across  this 
charge  distribution  then  aocounts  for  the  forward 
directivity  associated  with  the  relativlstlcally-movlng 
elec cron.  Thus  the  relativistic  nature  of  the  primary 
Compton  ourrents  oan  he  incorporated  directly,  through 
the  use  of  the  Lienard-Wieohert  fields.-  or  indirectly, 
through  the  careful  allowance  for  retardation  across  the 
full  EMP  source-current  distribution. 

The  treatment  of  the  secondary  eleotrons  comprises 
the  third  aspect  of  the  EMP  problem  which  has  been 
considered  in  the  previous  reports  in  this  series. 

These  secondary  eleotrons  are  produced  by  ionizations 
along  the  Compton  tracks,  and  by  other  ionization 
prooesses  such  as  the  photoelectrio  ejection  of 
eleotrons  from  air  molecules  by  incident  ultraviolet  * 
visible,  and  X-ray  photons.  The  secondary  eleotrons, 
once  released,  will  move  in  response  to  the  local 
electromagnetic  field,  and  will  in  turn  make  their 
contribution  to  this  field  distribution. 
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The  field  contributions,  made  by  the  secondary 
electrons  in  their  movement  in  the  local  fields, 
will  be  retarded-field  contributions,  Just  as  the 
fields  generated  by  the  Compton  electrons  were  retarded 
fields.  An  integral  equation  results,  whose  solution 
gives  the  distribution  of  the  secondary-electron  currents. 
Before  this  integral  equation  can  be  written  do*»m  in 
detail,  it  is  neoessary  that  expressions  be  given  for 
the  retarded  fields  associated  with  a  particular 
source-ourrent  distribution.  Once  these  expressions 
have  been  found,  the  integral-equation  problem  can  be 
given  explicitly,  and  solved  by  numerical  or  analytical 
methods. 

There  is  a  hidden  hazard  in  attempts  to  solve 
Maxwell* s  equations  by  numerical  methods  which  do  not 
incorporate  retardation  explicitly.  The  hazard  is 
attributable  to  the  fact  that  Maxwell *  s  equations  admit 
an  advanced-field  solution  in  addition  to  the  physically 
acceptable  retarded-field  solution,  A  general  solution, 
therefore,  will  be  a  superposition  of  retarded  and  advanced 
solutions,  unless  special  precautions  are  taken  to  exclude 
the  advanced  solutions  from  the  beginning.  As  will  be 
shown  later,  one  numerical  method  whloh  has  actually  been 
programmed  at  great  expense  violates  this  basic  physical 
causality  requirement.  Fortunately,  much  of  the  program 
oan  be  retained  when  the  non-causal  solution  is  converted 
to  a  causal  solution. 
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2.  GROWTH  OP  DIPOLE  MOMENT 

As  determined  from  two  successive  integrations 

o 

of  an  experimentally-measured  waveform,  the  dipole 
moment  associated  with  the  EKP  source- current  distribution 
has  a  time  dependence  which  shows  a  gradual  growth,  on 
a  time  scale  of  the  order  of  100  microseconds.  (In  some 
oases  there  is  an  overshoot,  but  this  also  is  gradual  on 
the  same  time  scale,  and  the  attainment  of . the  final 
dipole  moment,  after  the  moderate  overshoot,  coincides 
with  the  reduction  of  the  vertical  current  flow  to  zero.) 
This  100-nlcrosecond  time  scale  can  be  compared  with  the 
one-microsecond  time  scale  that  characterizes  the  actual 
nuclear  detonation.^  It  is  apparent  that  the  100-mlcrosecond 
radiated  waveform,  though  initiated  by  the  one-microsecond 
detonation,  does  not  have  its  time  dependence  determined 
by  the  time  scale  of  the  detonation's  chain  reaction,  but 
by  some  other  phenomenon  or  phenomena. 

After  the  prompt  effects  of  the  nuclear  chain  reaction, 

there  are  delayed  reactions,  in  which  gamma-rays  are  emitted 

from  the  detonation  products.  However,  these  delayed 

emissions  decay  with  a  time  constant  of  the  order  of 

relatively 

one  microsecond,  and  will  be/insignif leant  long  before  the 
lapse  of  100  microseconds. 
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The  prompt  Compton  currents,  associated  with  the 
prompt  gamma-rays  from  the  nuclear  chain  reaction,  will 
be  confined  in  time  to  an  interval  associated  with  the 
detonation  itself.  However,  the  propagation  time  must  be 
added,  since  the  Compton  ejections  at  a  distance  r  from 
the  detonation  center  will  occur  at  a  time  which  Is  later, 
by  the  time  interval  r/c,  than  the  moment  at  which  the 
ejecting  gamma-rays  were  emitted  by  the  nuclear  reaction 
at  r  »  0.  The  attenuation  in  air  of  the  pertinent 
gamma-rays  leads  to  an  attenuation  with  time  which  has 
about  the  same  time  constant,  one  microsecond,  as  the 
delayed  gamma-ray  emissions.  Thus  the  gamma-ray  propagation 
time  will  be  inadequate  to  account  for  the  slow  growth  of 
the  dipole  moment  associated  with  the  electromagnetic 
effects  of  a  nuclear  detonation  which  is  in  the  lower 
atmosphere,  near  the  ground.  . 

It  can  be  noted,  however,  that  the  range  of  the 
ultraviolet  radiation  from  the  detonation  is  substantially 
greater  than  the  range  of  the  high-energy  gamma-rays 
responsible  for  the  Compton  ejections.  Furthermore, 
the  forward  directivity  of  the  ejected  Compton  electrons 
leads  to  a  voltage  pulse  which  is  projected  outward  at  the 
velocity  of  light,  far  beyond  the  the  distance  reached  by 
the  Compton  electrons  themselves.  This  voltage  pulse 
has  the  polarity  to  drive  secondary  electrons  In  the  same 
radial  direction  in  which  the  Compton  electrons  were  moving, 
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voltage 

and  in  this  way  the/pulse  has  the  effect  of  extending 
greatly  the  range  to  which  the  upward  electron  motion 
reaches.  E|y  extending  the  range  of  the  vertical  current, 
this  hybrid  phenomenon  (radial  electric  field  from 
inner  Compton  electrons,  multiplied  by  conductivity 
from  outer  photoionization  processes)  provides  a  mechanism 
which  can  explain  the  observed  time  scale  f^r  the  growth 
cf  the  dipole  moment  associated  with  the  EMP  souroe- 
ourrent  distribution. 

In  addition  there  may  be  higher-order  Interaction 
processes;  in  which  the  secondary-electron  motion 
at  early  times  contributes  fields  which  lead  to 
secondary-electron  motion  at  later  times.  The  computation 
of  these  complicated  relaxation  processes  must  be  left  to 
analytic  or  numerical  methods,  and. cannot  readily  be 
foreseen  by  way  of  general  principles.  For  the  computation 
to  have  physical  significance,  however,  it  is  essential 
that  retardation  be  incorporated  correctly,  so  that 
effects  will  not  precede  their  causes,  and  the  propagation 
of  electromagnetic  fields  will  not  be  faster  than  the 
velocity  of  light. 


Page  11 


3.  RELATIVISTIC  EFFECTS 

2  *3 

In  earlier  reports  *  two  source  models  were 
analyzed  relativistloally  with  the  use  of  the  field 
expressions  obtained  from  the  Lienard.-Wiechort  potentials 
One  of  these  models,  called  the  shell  modex,  was  carried 
to  the  point  of  waveform  calculation.  In  tr.is  model, 
a  spherical  shell  of  electrons  is  ejected  at  one  radial 
distance,  moves  outward  at  a  relativistic  velocity,  and 
then  is  deposited  at  a  larger  radial  distance.  That  is, 
the  electrons  composing  this  moving  shell  have  all  been 
ejected  simultaneously  from  an  inner  spherical  surface, 
all  move  radially  outward  at  the  same  velocity,  and  are 
simultaneously  deposited  on  an  outer  spherical  surface. 
The  inner  shell  is  initially  uncharged,  so  that  when 
the  negatively-charged  electrons  are  ejected,  they 
leave  behind  a  stationary  shell  of  positive  ions. 

Two  different  symmetries  were  considered.  One  was 
spherical  symmetry,  with  no  angular  dependence  of  the 
density  of  electrons  on  the  moving  shell.  The  other 
was  cosine  asymmetry,  with  the  electron  density-set 
proportional  to  the  cosine:- of  the  polar  angle  on  the 
spherical  surface,  (This  makes  the  lower  hemisphere 
positively  charged,  with  positrons,  or  with  the  Images 
of  the  upper  electrons  in  a  horizontal  ground  plane.) 
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For  the  case  of  spherical  symmetry,  the  results  of 

the  relativistic  calculations  were  identical  with  the 

results  obtainable  from  Gauss* s  law  of  the  electrostatic 

» 

field.  Of  the  six  components  of  the  electric  and  magnetic 
fields,  only  the  radial  component  of  the  electric  field 
does  not  vanish  by  symmetry.  Furthermore,  for  an 
observation  point  outside  the  outer  spherical  surface, 
the  magnitude  of  this  radial  component,  Er,  vanishes  at 
all  times,  when  the  contributions  of  the  moving  electrons, 
the  positive  ions  left  behind,  and  the  stopped  electrons 
that  have  reached  the  outer  spherical  shell,  are  all 
added  together.  These  separate  contributions,  however, 
do  not  vanish  individually. 

It  is  instructive  to  examine  the  role  played  by 
retardation  across  the  source  region.  Because  of  the  • 
finite  propagation  velocity  of  electromagnetic  effects, 
an  observer  cannot  be  sure  that  the  shell-model  source 
is  actually  spherically  symmetric  until  enough  time  has 
elapsed  to  permit  him  to  receive  signals  from  the 
charged  particles  on  the  far  side  of  the  source  region. 
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Figure  1  shows  an  early  stage  in  the  development 
of  the  shell  model.  What  is  shown  is  the  source 
distribution  contributing  to  the  field  component  Er 
at  a  particular  observation  point,  as  the  source  would 
appear  to  that  observer  with  retardation  taken  into 
account,’  At  the  early  moment  depicted  here,  the  observer 
is  receiving  the  fields  emitted  by  a  sector  of  moving 
electrons,  on  the  near  side  of  the  shell-model  source, 
and  the  electrostatic  fields  from  the  positive  ions 
which  they  left  behind  on  the  inner  spherical  surface. 

It  is  too  early  for  the  observer  to  receive  the  electro¬ 
static  fields  from  any  of  the  stationary  negative  ions 
that  will  be  formed  when  the  movins  electrons  come  to 
rest  on  the  outer  spherical  surface. 

Figure  2  shows  a  somewhat  later  stage.  The  nearer*  ... 
electrons  have  come  to  rest,  forming  a  sector  of  negative 
ions  vrhich  is  bounded  by  the  dashed  lines  in  the  figure, 

A  band  of  moving  electrons  is  also  visible,  and  the 
exposed  positive  ions  on  the  inner  spherical  surface 
now  cover  a  sector  which  is  greater  than  a  hemisphere. 
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Figure  a..  The  shell  model,  as  seen  from  the  observation 
point  at  an  early  instant,  before  any  of  the  stopped 
electrons  (negative  ions)  can  be  seen. 
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Figure  J>,  The  shell  modelj  as  seen  from  the  observation 
point  at  a  moment  somewhat  later  than  the  moment  shown 
in  Fig,  I.  Some  of  the  electrons  have  reached  the  outer 
shell,  forming  negative  ions,  and  more  than  half  of  the 
positive  ions  on  the  inner  shell  have  been  exposed. 
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Figure  3  is  a  still  later  stage.  The  full  inner 
shell  of  positive  ions  is  exposed,  and  most  of  the 
ejeoted  electrons  have  reached  the  outer  shell  and  have 
ootne  to  rest  as  negative  ions.  A  small  sector  of  moving 
electrons  remains  visible,  but  these  will  shortly  reach 
the  outer  shell  and  stop  there. 

It  should  be  emphasized  that  the  time  scale  for 
these  three  figures  represents  the  observer* s  time  scale, 
and  the  source  distributions  are  those  that  he  vrould  sense. 
The  appearance  of  asymmetry  is  solely  due  to  the  differences 
In  propagation  time  from  different  parts  of  the. source  to 
the  observation  point.  The  source  distribution  itself  is 
actually  spherically  symmetric  in  this  Instance, 

For  this  shell  model  the  separate  contributions  of 
the  positive  and  negative  ions  and  the  moving  electrons, 
including  the  fields  generated  by  the  ejection  and 
deceleration  processes,  can  all  be  evaluated  in  closed 
form.  The  moving  electrons  are  treated  reiativlstlcally. 
When  the  contributions  are  added  together,  the  resulting 

magnitude  for  E  at  an  observation  point  outside  the 

*  • 

source  region  is  found  to  be  exactly  zero  at  all  times. 


Figure  3.  The  shell  model  at  a  late  stage.  All  of  the 
positive  Ions  on  the  Inner  shell  have  been  exposed,  and 
most  of  the  electrons  have  reached  the  outer  shell  and 
stopped  there  to  form  negative  Ions,  Because  of  the 
retardation,  a  few  of  the  moving  electrons  on  the  far 
side  of  the  source  remain  visible  at  the  observation 
point. 
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If  spherioal  symmetry  is  relinquished,  then  there 
Is  no  longer  this  exact  cancellation  of  the  different 
contributions  to  the  radial  electrio  field,  Er,  at  an 
observation  point  outside  the  shell-model  source 
distribution.  Also,  there  is  no  longer  the  limitation 
of  the  generated  electromagnetic  field  components  to  the 
single  component  Ey.  In  particular,  if  a  weight  function 
is  introduced,  weighting  the  currents  by  the  factor 
cosine  0 ,  where  0  is  the  polar  angle  measured  down 
from  the  positive  z-axls  (assumed  to  be  vertically 
upward) ,  then  the  result  is  a  source  model  which  will 
be  called  the  shell  model  with  cosine  asymmetry. 

A  closely  related  model  is  the  opposed-hemisphere 
shell  model.  For  this  model,  the  source  currents  can  be 
considered  to  be  the  upper  half  of  the  shell  model  with 
spherioal  symmetry,  together  with  the  image  currents 
which  would  accompany  such  a  hemispherical  shell  source 
if  it  were  located  dlreotly  above  a  perfectly  conducting 
ground  plane.  Figure  4  shows  the  distribution  of 
positive  and  negative  ions  in  the  opposed-hemisphere 
shell  model,  after  the  electrons  ejected  from  the  inner 
shell  (and  their  positively-charged  images)  have  been 
deposited  on  the  outer  spherical  shell. 


Figur-1  The  opposed-hemisphere  shell  model,  after  the 
ejected  electrons  (and  their  positively-charged  images) 
have  traveled  from  the  ejection  points  on  the  inner 
spherical  shell  to  the  deposition  points  on  the  outer 
spherical  shell. 
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When  the  opposed-hemisphere  weight  function  is 

expressed  as  a  superposition  of  spherical-harmonic 

inodes,  involving  weight  factors  which  are  Legendre 

polynomials  P  (cos  ©) ,  tlien  the  leading  term  has .the 

weight  function  for  the  sheli  model  with  cosine 

asymmetry.  This  weight  function  is  P^fcos©),  which 

* 

Is  Just  equal  to  the  cosine  of  ©.  For  the  shell 
model  with  cosine  asymmetry,  Figure  5  shows  the 
ion  distribution  after  the  moving  electrons  have  come 
to  rest  on  the  outer  spherical  surface. 

When  the  shell  model  with  cosine  asymmetry  is 
examined  in  detail,  it  is  found  that  the  time  history 
of  the  fields  at  an  observation  point  outside  the  outer 
shell  can  be  Interpreted  with  the  aid  of  Figures  1-3. 
However,  the  weight  function,  cosine©,  enters  as  a 
factor  multiplying  the  distributions  of  positive  ions, 
moving  electrons,  and  negative  ions.  The  symmetry  about 
the  z-axls  remains  in  this  model,  but  the  loss  of 
spherical  symmetry  means  that  the*e  are  three  nonvanishing 
field  components;  Ey,  EQ,  and  •  The  other  three  field 
components,  E^,  Hy,  and  H^,  vanish  as  a  consequence  of 
the  symmetry  about  the  z-axis  and  the  restriction  of 
the  electrons  (and  their  images)  to  radial  motion  only. 
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Figure  The  shell  model  with  cosine  asymmetry,  after  the 

ejected  electrons  (and  cheir  positively-charged  images) 
have .traveled  from  the  ejection  points  on  the  inner  spherical 
shell  to  the  deposition  points  on  the  outer  spherical  shell. 
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As  in  the  ease  of  spherical  symmetry ,  the  field 
contributions  of  the  ions  and  electrons,  at  an  outside 
observation  point,  can  be  evaluated  in  closed  form.  The 
result  for  the  case  of  Ej,  has  been  given  explicitly 
in  Reference  3.  When  the  contributions  vrere  added 
together,  it  was  found  that  for  an  observation  point 
above  the  source  (on  the  positive  z-axis  at  a  position 
above  the  outer  shell),  the  radial  electric  field  rapidly 
approached  a  negative  peak,  then  more  slowly  decoyed  back 
to  a  residual  value  associated  with  the  static  d:-pole 
moment  of  the  charge  distribution  in  Figure  5*  The 
polarity  of  Er  was  at  all  times  negative,  in  a  direction 
to  drive  secondary  electrons  upward  if  there  had  been 
any  secondary  electrons  present,  in  the  region  above 
the  shell -mo del  source  distribution. 

For  an  observation  point  located  along  the  .z-axis, 
the  field  components  EQ  and  vanish  by  symmetry. 

However,  for  an”  observation' point  located  in  the 
neighborhood  of  the  ground  plane,  to  one  side  of  the  source 
E^  and  have  their  maximum  values  while  Ep  vanishes. 

For  an  intermediate  location,  as  .in  Figures  1-3*  *11 
three  components  are  nonvanlshing  at  the  observation  point. 
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While  the  electrio  field  components,  Er  and  EQ, 

are  found  to  have  residual  magnitudes  associated  with  the 

dipole  moment  of  the  final  ion  distribution  (Figure  5), 

the  magnetic  field  component,  H^,  is  found  to  be  directly 

associated  with  the  moving  electrons  that  are  indicated 

in  Figures  1-3  (as  modified  by  the  cosine  weight  function). 

The  magnetic  field  is  zero  at  the  observation  point 

up  until  the  moment  when  the  first  of  the  moving  electrons 

becomes  1 visiDle'  at  this  observation  point,  and  is  zero 

again  for  times  later  than  the  moment  at  which  the  farthest 

moving  electron  comes  to  rest  on  the  outer  shell,  as  sensed 

at  the  observation  point  with  retardation  appropriately 

incorporated.  These  moments  marking  the  beginning  and 

defined  as 

ending  of  the  magnetic-field  waveform  are/the  initial  moment 
slightly  preceding  the  moment  depicted  in  Figure  1,  and 
the  final  moment,  Just  after  the  moment  depicted  in 
Figure  3. 

For  times  which  lie  between  the  initial  and  final 
moments,  the  magnetic  field  at  the  observation  point  can 
be  described  as  the  sum  of  two  terms.  One  of  these  terms* 
the  induction  term,  falls  off  with  radial  distance  in 
proportion  to  r  .  The  other  term,  the  radiation  term, 
falls  off  with  radial  distance  in  proportion  to  r”1. 

Each  term  vanishes  separately  for  times  which  lie  outside 
the  interval  between  the  initial  and  final  moments,  defined 


above. 
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4.  RETARDATION  AND  CAUSALITY 

The  shell-model  calculations  described  in  the 
previous  Section  treated  the  moving  electrons  as 

relativistic  parti-clec,  and  utilized  the  Lionard-Wiechert 

4 

field  expressions.  The  same  waveform  expressions 
can  be  obtained  from  a  d.istributed-source-cv.rrent 
picture,  replacing  the  moving-charged- particle  picture, 
provided  that  the  current  distribution  is  inserted  into 
an  integral  formulation  of  the  fielos^  in  which 
retardation  is  accurately  Incorporated, 

When  this  integral  method  is  used,  the  shell-model 
current  distribution.,  when  reduced  from  four-dimensional., 
space-time  to  the  twe -dimensional  (r,t)-plane,  has  the 
form  of  an  elongated  delta-function,  as  illustrated  in 
Figure  6.  For  this  reduction,  the  angular  dependence 
of  the  fields  and  the  current  has  been  expressed  in  terms 
of  the  appropriate  vector  spherical  harmonics^,  and  only 
the  dependence  upon  radial  distance,  r,  and  time,  t, 
remains  in  the  reduced-  problem. 


Figure  36,  The  problem  domain,  in  the  (r,t) -plane,  for  the 
shell  model  with  cosine  asymmetry.  The  source  current  is 
an  elongated,  tilted  delta-function,  shown  here  as  an  elongated 
oval.  The  region  in  which  the  magnetic  field  differs  from 
zero  is  bounded  by  the  two  dashed  lines,  and  by  portions  of 
the  lines  r  =  0  and  r=  ct. 


Page  26 


The  choice  of  the  origin  of  time,  in  Figure  6»  has 
been  made  in  such  a  way  that  the  ejection  of  the  electrons 
from  the  inner  shell  could  have  been  through  Compton 
processes  involving  gamma  rays  which  left  the  point 
r  *  0  at  the  time  t  «  0,  Thus  the  ejections  take  place 
at  jl  point  in  the  diagram  which  lies  on  the  line  r«ct, 
and  the  current  flow  extends  to  the  right  of  this  line. 

The  slope  of  the  elongated  delta-function  is  equal  to 
the  ratio  v/c,  where  v  is  the  velocity  of  the  moving 
ele. tron  in  the  shel] -model  picture,  and  c  is  of  course 
the  velocity  of  light, 

Figi  6  shows  the  region  where  the  magnetic  field, 
Hjj,  differs  from  zero.  For  an  observation  point  at  a 
particular  radial  distance,  r,  outside  the  source  region, 
a  horizontal  section  of  Figure  6  shows  that  the  magnetic 
field  will  differ  from  zero  for  a  finite  tine  interval. 
This  is  the  Interval  during  which  moving  electrons  are 
'visible*  at  that  observation  point. 

In  the  integral  method^  the  field  components  at 
a  space-time  point  P  are  expressed  as  explicit  integrals 
over  the  causally  accessible  source -current  distribution. 
Through  separation  into  vector-spherical-harmonic  modes, 
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the  problem  Is  reduced  to  two  dimensions,  and  the 
« 

field  components  for  a  particular  mode  are  expressed 
as  integrals  in  the  (r,t) -plane,  over  the  current 
components  for  the  same  mode. 

The  realm  of  this  integration  in  the  (r,t) -plane 
is  shown  in  Figure  7«  It  is  assumed  that  no  source 
currents  flow  in  the  region  to  the  left  of  the  diagonal 
line  r»ct,  because  it  is  the  EMP  problem  that  is 
to  be  treated,  and  the  nuclear  detonation  initiates  the 
source-current  flow.  The  currents  which  contribute  to 
the  field  components  at  a  space-time  point  P  are  then 
the  currents  in  the  shaded  regions  in  Figure  7.  The 
horizontal  shading  indicates  the  currents  which  can 
make  inductive  and  radiative  contributions  to  the  field 
components  at  P.  (Thus  the  magnetic  field  components 
at  P  must  be  generated  by  currents  in  the  region  that 
is  horizontally  shaded,)  In  addition,  there  are 
electrostatic  contributions  by  the  currents  in  the  rectangular 
region  with  horizontal  shading,  and  also  by  the  currents 
in  the  triangular  region  with  vertical  shading. 

It  is  the  limitation  Imposed  by  the  causality  diagram 

the 

of  Figure  7,  which  accounts  for/limited  region  of  nonzero 
magnetic  field  shown  in  Figure  6,  from  the  point  of  view 
of  the  integral  method. 


Figure  17.:,  The  causality  diagram  for  the  integral  method. 

The  field  components  at  the  space-time  point  P  are  determined 
by  the  currents  in  the  shaded  region.  The  inductive  and 
radiative  contributions  are  made  by  the  currents  in  the 
rectangular  region  shown  with  horizontal  shading.  The 
triangular  region  with  vertical  shading  produces! only  r 
electrostatic-field  contributions. 
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5. . SECONDARY  ELECTRONS 

The  treatment  of  the  secondary  electrons  comprises 
the  third  of  the  three  aspects  of  the  EMP  problem  that 
were  referred  to  in  the  Introduction.  Some  of  the  secondary 
electrons  are  those  that  are  released  by  the  ionizing 
collisions  that  serve  to  slow  the  Compton  electrons 
in  their  passage  through  the  air.  Other  secondary 
electrons  are  released  from  air  molecules  through  the 
photoelectric  ionization  processes  in  which  the  energy 
for  ionization  is  provided  by  X-rays,  ultraviolet  light, 
and  visible  light  from  the  nuclear  detonation.  The 
secondary  electrons  released  by  ionizations  along  the 
Compton  tracks  will  be  mainly  confined  to  the  inner 
regions  where  there  are  many  Compton- electron-producing 
gamma  rays.  The  X-ray  ionization  will  also  have  a 
relatively  short  range.  However,  the  ionization  which 
is  produced  by  ultraviolet  light  will  be  spread  over  a 
large  volume,  because  of  the  relatively  long  range  of 
ultraviolet  photons  in  air,  as  compared  with  X-rays  and 


gamma  rays. 
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The  secondary  electrons ,  once  released  by  whatever 
ionizing  process,  will  move  in  response  to  the  local 
electric  field  strength.  Their  motion,  in  turn,  will 
constitute  a  current  which  will  generate  its  own  electric 
and  magnetio  fields.  Since  all  of  the  processes  arc 
transient1,  and  since  there  are  propagation  delays 
involved  each  time  a  moving  electron  affects  another 
electron,  it  is  possible  in  principle  to  set  up  a 
mathematical  iterative  procedure,  beginning  at  very 
early  times  when  the  secondary  ionization  is  small. 

For  the  initial  steps  of  the  iteration  procedure, 
what  matters  is  that  the  product  of  the  conductivity 
associated  with  the  secondary  ionization,  and  the  electric 
field  associated  with  the  primary  Compton  current, 
should  give  a  secondary  current  which  is  small  in 
comparison  with  the  primary  Compton  current.  The  iteration? 
then  procoeds  by  steps  in  the  (r,t) -plane  which  are 
small  enough  so  that  numerical  instabilities  are 
avoided  and  a  smooth,  gradual  function  for  the  net 
ourrent  (the  sum  of  the  Compton  current  and  the  secondary 
Current).  Is  obtained. 
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It. is  also  possible  to  call  upon  analytical  methods 
for  the  solution  of  the  integral  equation,  which  gives 
the  motion  of  the  secondary  electrons  and  the  fields 
they  generate, ?  in  regions  of  space- time  where  the 
secondary  conductivity  is  not  very  great,  an  analytical 
iteration  procedure  can  be  used,  where  the  secondary 
current  is  written  as  a  sum  of  terms.  The  first  term 
is  the  product  of  the  conductivity  function  and  the 
electric  field  generated  by  the  Compton  electr^"  alone. 

The'  second  term  is  the  product  of  the  conductivity  function 
and  the  electric  field  generated  by  the  first  term  alono, 
and  so  forth.  This  approach  should  be  useful  for  early 
times  and  for  large  radial  distances. 

In  regions  of  space-time  where  the  secondary  conductivity 
is  very  high,  the  above  iteration  procedure  can  lead  to 
numerical  oscillation  or  divergence.  Here  the  straightforward 
numerical  approach  can  be  used,  but  an  alternative  is  one 
in  which  the  net  current  or  the  net  local  electric  field 
is  expanded  in  a  finite  set  of  linearly  independent  functions 
which  span  the  region  of  interest  in  the  (r,t) -plane. 

The  integral  equation  can  then  be. made  to  yield  a 
least-squares  condition  on  the  coefficients  in  the  above 
expansion. 
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6,  CONCLUSIONS 

Special  emphasis  has  been  laid  on  the  causality- 

requirement  ,  in  the  solution  of  the  EMP  problem.  The 

reason  for  this  emphasis  is  that  the  source  volume  is 

the  time 

large  enough  to  make/retardation/across  this  volume 
a  significant  part  of  the  generated  electromagnetic 
signal  waveform . duration,  A  solution  to  Maxwell's 
equations  which  does  not  ailow  properly  for  this 
retardation  may  have  little  or  no  relationship  tn 
the  physical  processes  which  take  place  in  the  vicinity 
of  a  nuclear  detonation. 

It  was  noted,  for  example,  that  the  shell-model 
calculation  showed  an  early-time  radial  electric  field 
which  had  the  polarity  to  drive  secondary  electrons 
upward  above  the  detonation  center,  and  that  this 
early-time  radial  field  coincided  in  time  with  the 
long-range  ionization  by  ultraviolet  light  from  the 
detonation,  and  was  therefore  important  in  the  analysis 
of  the  mechanisms  contributing  to  the  establishment  of 
the  large  dipole  moment  which  regains  after  the  EMP 
waveform  is  over. 
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There  Is  a  non- causal  numerical  solution  of 
Maxwell's  equations  which  has  been  applied  to  this 

EMP  problem.^  As  evidence  that  this  non-causal  solution 

8 

does  indeed  give  wrong  answers,  it  has  been  noted  that 
the  early-tine  radial  electric  field,  obtained  from 
this,  non-causal  solution,  has  the  wrong  polarity.  Instead 
of  driving  secondary  electrons  upward  above  the  detonation 
center,  the  early-time  radial  electric  field  from  the 
non-causal  calculation  drives  secondary  electrons  downward. 
Thus  an  important  physical  process  contributing  to  the 
dipole-moment  establishment  does  not  appear  in  this 
non-causal  solution. 

It  is  fortunate  that  the  non-causal  iteration  method 
can  be  separated  from  the  remainder  of  this  numerical 
solution  of  Maxwell's  equations,  and  can  be  replaced  by 
a  causal  iteration  procedure,  without  changing  the  portions 
of  the  program  which  deal  with  electron  attachment  and 
air  chemistry,  and  v’  th  the  gamma  rays  and  the  primary 
source  current. 
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|  ABSTRACT 

In  analogy  with  the  familiar  addition  formula  for 

» 

Legendre  polynomials,  a  generalized  addition  theorem  is 
proved,  A  general  spherical  harmonic,  depending  on  the  two 
angles  9^  and  ,  is  expressed  as  an  expansion  Involving 
spherical -harmonic  functions  of  (©,p)  and  of ‘ (©'  ip* )•  Die 
six  angles  are  related  to  each  other  through  the  equations 
cos  ©^  =  cos  9  cos  ©'  -  sin  9  sin  9*  cos  (p*  -  p  cos  9)  , 
sin  cos(p1-p)  =  sin©  cos  9*  +  cos  ©sin  ©*  cos  (p*- p  cos  9)  , 
sin  ©^  slnCp^-p)  *  sin©'  sin(p*-p  cos  9)  ,  The  expansion 
is  then  used  in  the  proof  of  an  integral  theorem,  for 
spherical  harmonics. 


j 

\ 
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1.  INTRODUCTION 

The  familiar  addition  formula  for  Legendre  polynomials* 
can  bo  written  as 

Pn(cos6f)  =  Fn(cos  ex)  Pn(cos  0)  ■' 
jx-n 

+  .2^  ~  -  p£(cos  ©-,)  p£(cos-0)  cos^f^-#)  ,  (1.1) 

t— »  (r+m)I 
1 

where  •  * 

cose1  s=  cos  cos  6  +  sin  6^  sin  0  cos  (0^-0)  •  (1.2) 

The  angle  0*  can  be  interpreted  as  the  included  angle 
between  the  two  vectors  rn  and  r,  whose  directions 

e»J,  «** 

are  specified  by  (6^*^)  and  (6,0),  respectively. 

The  Legendre  polynomials  are  a  subset  of  the 
more  general  spherical  harmonics.  Eq.  (1.1)  is  a 
relationship  between  members  of  such  a  subset,  on  the 
left-hand  side,  and  members  of  two  full  sets,  on  the 
right-hand  side.-  It  is  to  be  expected  that  other 
similar  relationships  exist,  in  which  a  full  set  of 
spherical  harmonics,  Involving  0*  and  a  suitably 
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two 

defined  0X ,  is  related  to  the/full  sets  of  spherical 
harmonics  involving  ( , ^ )  and  (6,0). 

It  is  first  necessary  to  find  a  suitable  definition 
for  the  angle  0'  which  is  to  accompany  ©*.  A  simple 
definition  is  found,  and  the  reason  for  its  selection 
is  explained.  An  equation  of  the  form  (1.1)  is  then  given, 
in  which  P^(cos  6^)  is  expressed  a3  a  summation  over  spherical 
harmonics  which  are  functions  of  0,0,6* ,0' .  Finally,  this 
equation  is  generalized  so  that  a  general  spherical  harmonic, 
written  as  a  function  of  is  then  given  as  a  summation 

over  spherical-harmonic  functions  involving  6,0,6' ,0' . 
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2.  DEFINITION  OF  0* 

The  vectors  r  and  r^,  with  polar-coordinate 
components  (r,e,$)  and  are  illustrated 

in  Figure  1.  These  vectors  are:  shown  with  respect  to 
a  rectangular  coordinate  system,  with  axes  labeled 
by  x,  £,  and  z. 

Figure  2  shows  a  new  set  of  coordinate  axes, 
labeled  x' ,  y  * ,  and  z' ,  which  are  obtained  from  the  first 
set  by  two  rotations.  The  first  is  a  rotation  about  the 
z-axis  by  the  angle  p;  this  rotation  moves  the  y-dlrection  into 
the  position  of  the  y'-axis,  which  thus  must  lie  in  the 
(x,y) -plane.  The  second  rotation  is  about  the  .  '  y'-axis, 

and  moves  the  z-direction  down  by  the  angle  0,  until  it  lies 
in  the  position  of  the  z'-axis,  which  lies  along  the  vector  r. 

In  this  way,  by  the  two  rotations  through  the  angles  and  0, 
the  (x,y,z) -directions  are  moved  into  the  positions  shown 
as  the  (x* ,y* ,z' )-axes. 

Figure  2  thus  defines  the  z'-axis  unequivocally,  since 
this  axis  must  lie  along  the  direction  of  the  vector  r.  The 
y'-axis ’will  also  be  defined  uniquely,  if  ©  is  greater  than 
zero  and  less  than  tr.  The  plane  through  the  origin  of 
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to  the  (x,7,z)-axes. 


shown  in  relation 
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coordinates  which  Is  perpendicular  to  r  will  then  intersect 
the  (x,y)-plane  In  a  straight  line  which  contains  the 
y'-axis.  The  direction  of  this  axis  can  then  be  obtained, 
as  shown  in  Fig.  2,  through  the  application  of  the  right-hand 
rule.  If  a2  is  a  unit  vector  In  the  direction  of  the  z-axis, 
while  r/r  is  a  unit  vector  in  the  direction  of  the  z’-axis, 
then  the  unit  vector  in  the  y’-directlon  is 


»zxr 

r  sin  Q 


(2.1) 


If  the  vector  r  happens  to  be  parallel  to  the  z-axis, 

.  ao  that  sin  Q  is  equal  to  zero  and  the  vector  product 
a^xr  also  equal  to  zero,  then  the  direction  of  the 
y'-axis  becomes  indeterminate.  Since  any  value  of  can 
he  used  in  specifying  the  direction  of  r,  when  €  is  equal 
to  zero  or  to  tt,  It  Is  apparent  from  Fig.  2  that  the 
y'-axis  can  be  directed  in  any  direction  within  the  (x,y)-plane. 

•It  Is  essential  that  provision  be  made  for  this 
indeterminacy  In  the  orientation  of  the  y*-axls  (and  the 
tomes  ponding  indeterminacy  in  the  orientation  of  the 
x'-axls)  when  the  angle  is  defined,  to  accompany  the 

angle  ©' . 
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What  is  desired  is.  that  the  definition  of  0X  should 
permit  the  specification  of  the  direction  of  the  vector  r, 

— i. 

in  terms  of  the  direction  of  the  vector  r  and  the  direction 

09 

angles  .  ©'  and  0y •  In  the  ordinary  situation,  in  which  • 
r  is  not  parallel  or  antiparallel  to  the  z-direction, 

09 

the  y'-axis  is  uniquely  defined  and  no  problem  arises. 

However,  when  r  lies  in  the  positive  or  negative  z-direction, 
the  angle  0  is  not  well  defined  and  it  becomes  difficult  to 
give  the  azimuth  angle  0^  for  the  vector  r^,  in  terms  of  the 
ill-defined  azimuth  angle  0  of  the  vector  r. 

09 

The  solution  of  this  problem  of  azimuth  Indeterminacy 
is  shown  in  Fig.  3.  With  respect  to  the  (xf ,y ' ,z ' )-axes, 
the  direction  of  the  vector  is  specified  by  the  polar  ' 

angle  ©*'  and  by  an  azimuth  angle  which  is  given  by  the 
expression  (0*  -0  cos  Q) »  When  r  is  not  parallel  (or 
antiparallel)  to  the  z-axls,  then  the- use  of  (0*  -0cos©) 
Instead  of  0 1  represents  a  simple  displacement;.of :■  the  . 
azimuth-angle  coordinate.  However,  in  the  exceptional 
cases,  where  r  is  parallel  [or  antiparallel!  to  the 
z-axls,  the  r  ;  ichoice, ‘  .  of  (0*  -0cos©)  provides  an 
unequivocal  specification  of  the  azimuth  angle  0^,  even 
when  0  Itself,  and  0 1  accordingly,  are  both  indeterminate. 


i 
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3-  UK IT  VECTORS 

In  terms  of  the  rectangular  set  of  unit  vectors, 

1 ]  J,  k,  directed  'along  the  (x»y»z)-axea,  there  are 

mm  mm  mm 

three  polar- coordinate  unit  vectors  associated  with  the 
vector  r  ;  these  are : 

a  =  1  sin  ©cos  0  +  J  sin  ©  sin#  +  kco3  ©  t 

a©  =  i  c°s©  cos  0  +  J  cos  ©sin  /  -  kain©,  (3*1) 

aw  *  -  lain#  +  Jcos#. 

mo  mm 

A  similar  set  of  unit  vectors  is  associated  with  the 
vector  r,  : 

a^r  *  1  sin  ©^^  oosp^  +  J  sin  8^^  sin  +  k  cos  ©^  , 

he  B  1  cos  ©1  cos  pl  +  J  cos  8.^  sir*  pj^  -  ksin©^,  (3*2) 

Sip  *  -  i  sinp1  +  icosPj^, 

The  primed  axes  In  Pig,  2  have  been  chosen  to  lie  in  the 
directions  of  the  unit  vectors  (3.1) »  so  that  the  primed 
unit  vectors,  i* ,  J1,  k' ,  are  given  by 
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Ihe  unit  vector  which  is  parallel  to  the  vector  r^ 

©aii  be  written  as  alr  ,in  the  form  given  in  (3.2),  but  it 
c&ii  .s3Uis  be  expressed  in  terms  of  the  primed  axes; 

a^r  **  i*  sin  ©' cos(p' -P  cos  ©)  +  J '  sin  ©'  sin(p'  -p  cos  ©) 

+  k'  cos  ©'  *  (3.^) 

mm 

Equations  (3.1-*0  yield  the  trigonometric  transformations. 

cos  6^  *  cosQcos©*-  sin&sin©'  cos(p'-p cos  ©)  ,  (3.5a) 

•in  ©^  cos  (p^-p)  *  sin©oos©*  +  oos©sin©'  cos(p*  -pcos©),  (3.5to) 
•in  ©1sin(p1-p)  =  sin  ©•  sln(p* -p cos  ©)  ,  (3.5®) 

and  the  inverse  transformations 

cos  ©'  «  cos  ©1  cos  ©  +  sin  ©1  sin  ©  cos  (p^-p)  ,  (3.6a) 

•in©'  cos (p* - p  cos  ©)  »  -  cos  ©^sln©  +  sin  ©^  cos  ©  cos  (p^-p)  ,  (3.6b) 

sine’  slntp'-p  cos  6)  =  sin  9,  sln(p^-p)  , 


(3.6c) 
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4.  ADDITION  THEOREM 

The  addition  formula  (*1.1)  is  based  on  (3.6a),  but 
a  similar  formula  can  readily  be  written  down,  based  on 
(3.5a).  t 

Pn(cos  *  Pn(cos  0)  Pn(cos  6* ) 

*  +  2  )  p^(cos  ©)  P^(cos  ©• )  (-1)M  cosf/^p'-p  cos  6)1.  (4.1) 

(n+/t) !  n  n  L  J 

^=1 

It  is  this  formula  which  will  be  generalized  into  an 
addition  theorem  applicable  to  any  spherical  harmonic, 
not  Just  to  the  Legendre  polynomial  Pn(cos  ©1)  . 

2 

In  the  notation  of  Morse  and  Feshbach  ,  a  general 
spherical  harmonic  ban •' be  ‘  defined  by 

=  exp ( lnpx 5  pN(  cos©1)  ,  (4.2) 

where 

riPkcose,)  =  (sine,)1”1-— i-.p_(eose1) .  (4.3) 

"  1  1  4(cos  e^1”1  n  A 

It  will  be  assumed  that  0,  p,  6* ,  p%  are  independent 
variables,  while  and  p^  are  the  dependent  functions 
given  by  (3.5).  Two  operators,  and  D~  ^  ,  will 
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be  defined  by 

T  - 1  a  pa-  i 

®»P  L  s*n  ®  a®  oos  ^  ®P  sinz9  cos  9 

Ll  terms  of  these  operators,  the  spherical  harmonic  (4,2) 
can  be  written  as 


-  eip(lmp)  sltiNe  (De>)5)m  ^(cose.^,  (4.5) 


m 


in  which  the  operator  expression  (Da  »)  has  the  meaning 
^D9,p'  when  m  is  positive,  (D^p)^  when  m  is  negative. 
The  proof  of  (4,5)  is  based  on  the  easily  verified  results, 


.  exp(iCn)  sin  9, 

DJ  0  cos  9,  *  - i - *■ 

1  exp  (ip)  sin  9 

(D^p^oos  0-l  *  0, 


(4>6a) 


(4.6b) 


and  the  complex  conjugate  equations  involving  • 

It  is  apparent  also  that 

)£(e,p)  -  eip(inp)  sln,m|9  (DeijS)“Pn(cose).  (4.?) 


When  Pa(cos  9^)  bn  the  right-hand  side  of  (4,5) 
Is  replaced  by  the  series  expression  (4,1) ,  the  result 
is  the  generalized  addition  theorems 
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*  x£(®.p)  Pn(cos6')  +  2^ 


*=n 

+  2  y  fe£2i  (.!,»*  p^oose-) 


/C*l 


(n+>i)! 


•  exp(imp)  sin^e  (De^)n^(oos  e)  cos^>i(p'- p  cos  9)]}.  (4.8) 

Equation  (4.8),  together  with  the  transformation  equations  (3*5)» 
constitutes  the  desired  generalization  of  (1*1)  and  (1.2)* 

The  use  of  the  azimuth-angle  expressions  (p^-p)  and 
<P*-P  cos  0)  is  a  necessary  complication,  without  which  the 
important  equation  (4.5)  could  not  be  established.  However, 
once  a  choice  of  (n,m)  has  been  made,  and  the  operations 
in  (4.8)  carried  out,  the  angle  (p*-p  cos  6)  can  be  replaced 
by  a  re-defined  p' ,  here  and  in  Fig.  3,  if  this  is  desired 
in  a  particular  application  of  the  formula.  The  replacement 
cannot  be  introduced  until  after  the  operations  (D^p)”1  have 
been  completed. 
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^  5.  INTEGRAL  THEOREM 

The  addition  tneorem  (4.8)  leads  directly  to  an 
integral  theorem,  Eq,.  (4.8)  expresses  the  spherical 
harmonic  as  a  function  of  four  independent 

variables,  ©,p,©* ,p* .  The  dependence  upon  p*  is 
particularly  simple,  as  can  be  seen  when  (4.8)  is 
written  in  the  form 


/l  =n 


3^(0!, Pj.)  =  x“(e,p)  Pn(eos  e')  +^[ac  eosQip*)  +  Ag  sln(|ip')] 

_P-*=1 


(5.1) 


where  the  coefficients  AQ  and  A  depend  upon  the  variables 
e,p,e<  and  the  parameters  n,  m,/j,  but  not  upon  p* . 

Equation  (5*1)  can  now  be  integrated  with  respect  to  p* , 
with  the  other  three  variables  held  constant.  Each  term  In 
the  summation  over  yji  integrates  to  zero,  leaving  only: 


(®1Pp^)  dp< 


-  2ir  ^"(e.p)  Pn(oos  e*j . 


Other  related  integral  theorems  are  given  in  a  separate 
article^. 


(5.2) 
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ABSTRACT 

With  the  aid  of  a  generalized  addition  theorem 
for  spherical  harmonics,  previously  obtained,  six 
integral  theorems  for  vector  spherical  harmonics  are  proved. 

A  source-point  direction,  is  first  expressed  in 

terms  of  a  field-point  direction,  (0,p),  and  a  polar-coordinate 

angle-pair ,  (©* ) ,  which  has  as  its  polar  axis  the 

field-point  direction.  For* a  particular  choice  of  n'and  m,  all  ths 

components  of  the  vector  spherical  harmonics  for  the  source, 

expressed  in  terms  of  (0-^,^),  are  integrated  over  the 

relative  azimuth  angle,  p1 , while  the  field-point  direction, 

(©,/>),  and  the  relative  polar  angle,  ©* ,  are  held  fixed.  .•  The 
result  in  each  case  is  a  spherical  harmonic  or  vector  spherical 
harmonio  of  the  field-point  direction,  with  the  same  n  and  m 
but  now  depending  on  (©,£>)  instead  of  (0^,^),  multiplied 
.by  an  explicit  function  of  the  relative  polar  angle,  ©• • 
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1.  INTRODUCTION 


In  the  preceding  article'*' ,  a  generalized  addition 
theorem  for  spherical  harmonics  was  established.  This 
addition  theorem  led  directly  to  the  proof  of  an  integral 
theorem  for  spherical  harmonics.  In  the  integration, 
a  vector  is  swung  about  a  fixed  vector  r,  with  a 

constant  angle  ©’  separating  the  two  vectors.  A  suitably 
defined  azimuth  angle,  locates  the  azimuth  of  In 

its  motion  about  r.  This  integration,  "over  a  full  circuit, 
gave  the  result: 


*S'=ar 

I  =  2n  xj(e,(()  pn(cos  e') , 


(l.i) 


j!'“o 

vhere  ^  U  .  general  spherical  harmonic,  and  Pn  .".f; 
is  a  Legendre  polynomial. 


In  this  article  the  result  (1*1)  will  be  generalized 
to  apply  to  problems  which  arise  in  the  use  and  application 
of  vector  spherical  harmonics,  vector  functions  which  are 
closely  related  to  the  scalar  spherical  harmonics  . 
Because  there  are  three  independent  vector  spherical 
harmonics  for  each,  choice  of  (n,m),  and  each  of  these  has 
three  vector  components,  the  generalization  of  (1.1)  leads 


to  nine  Independent  scalar  equations.  However,  these 
nine  scalar  equations  group  naturally  into  three 
radial-component  equations,,  each  of  scalar  form, 
and  three  transverse- component  equations,  each  of  which 
combines  two  scalar  equations  into  a  vector  equation  with 
only  two  independent  components  rather  than  three. 

There  are  thus  only  six  separate  integral  theorems 
to  be  considered,  rather  than  nine.  This  article  gives 
the  proof  of  each  of  these  six  integral  theorems. 
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2.  VECTOR  SPHERICAL  HARMONICS 


The  notation  for  the  vector  spherical  harmonics, 

to  be  used  in  this  article,  follows  that  used  by  Morse 
o 

and  Feshbach  .  The  scalar  spherical  harmonics  are  defined  by 


j£(0.p)  =  ezp(imp)  p£“'(cos  ©)  , 


(2.1) 


where 


'oos  ©) 


(sinO),ml - ~ - rr,  Pn(cos  ©)  . 

d(oos  0),nl  n 


(2.2) 


The  vector  spherical  harmonics  separate  into  three  sets; 


P“(6,P)  =  ar  3^(0, J3)  , 


B”(©.P)  - 


Mn+iifr 

+  1)  sine  L 


(2.3) 


(n-tmf+1)  m  ^  (njjmj)  m 

(n+1)  n  a-l_ 


Ym 


1  m  (2n+l) 
+  ~P  n(n*M) 


1 


(2.^) 


C“(9,P)  = 


Vn(n+i) 
(2n+l)  sin  6 


|se  “ 


m  (2n-H)  m 
n(n+l)  ^ 


-  AP[ 


(n-M+1)  m  (n+|mp 

(n+1)  n+1  "*•  n 


•  +  ND  yQ  ll 

“  JJ 


(2.5) 


In  (2,4)  and  (2,5)  the  arguments  of  the  scalar  harmonics  are 
(6,p)  in  each  case.  Explicit  forms  for  the  unit  vectors  ap  , 
a^  ,  and  a.p  are  given  In  Eqs.  (3.1)  of  Reference  1. 
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2.  VECTOR  SPHERICAL  HARMONICS 

The  notation  for  the  vector  spherical  harmonics, 

to  be  used  in  this  article,  follows  that  used  by  Morse 
o 

and  Feshbach  .  The  scalar  spherical  harmonics  are  defined  by 
j£(e,p)  =  erp(imp)  P^kcos©)  ,  (2.1) 

where 

|j*j!  — ' 

pN(oose)  =  (sin«),m| Pn(oos6)  .  (2.2) 


The  vector  spherical  harmonics  separate  into  three  sets: 


P*  C©,p)  =  ®r  (G’?>  * 


Vn(n +TY  f 
(2n  +1)  sin 0  . 


(n-Imf+I)  m 
(n+1)  .  n+l 


(2.3) 


(n+]m|)  vm 


n 


+  a 


C“(9,p)  = 


Vn(n+1 


(2n+l)  sin 


) _  C  1  a  (2n4 

n©|~e  n(n+l) 


1  m  (2n+l) 
n(n*HL) 

(2n+l)  m 


1 


V-i 


(2.*0 


~*p[ 


(n~|mHl)  m  (gjN)  Ym 

(n+1 )  **+1  ’•  n  **-1 


(2.5) 


In  (2,4)  and  (2,5)  the  arguments  of  the  scalar  harmonics  are 
(©,p)  in  each  case.  Explicit  forms  for  the  unit  vectors  ap  , 
,  and  a ,p  are  given  in  Eqs,  (3.1)  of  Reference  1. 
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*  .  3.  THEOREM  I 

The  s lx  integral  theorems  involve  integrations 
over  p*  of  components  of  the  vect spherical  harmonics 
which  are  expressed  as  functions  of  For  example, 

the  vector  spherical  harmonic  JP°  (6^  ,p^)  has  the  form 

Jn<el'V  -  8ir£<W*  O.D 

obtained  from  (2,3),  where  the  unit  vector  a^r  ,  and  the 
orthogonal  unit  vector?  and  a^  ,  have  been  given 

explicitly  in  Eqs,  (3.2)  of  Reference  1.  The  geometrical 
configuration  of  the  vectors  r  and  r^  is  shown  in 
Figs,  1-3  of  Ref.  1,  while  the  trigonometric  transformations 
are  given  in  Eqs,  (3.5)  and  (3.6)  of  Ref,  1, 

In  the  integration  over  p* ,  the  vector  r  is  held 
fixed,  while  the  vector  r^  is  swung  In  an  arc  with  the 
angle  ©*  held  constant.  The  path  of  integration  is  shown 
here  in  Fig.  1.  Of  the  four  independent  angular  variables, 
e,p,0* ,p’ ,  only  p*  varies,  but  both  of  the  angles  ©1 
and  p^  will  vary  during  the  integration.  These  latter 
are  treated  as  dependent  variables,  through  the  use  of 
Eqs.  (3.5)  of  Ref.  1, 

The  direction  of  the  vector  function  (3.1)  changes 
during  the  integration  over  p* ,  since  a-^r  is  a  function  of  p 
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The  vector  function  (3.1)  will  therefore  be  resolved  Into 
Its  components  before  the  Integration,  and  the  simplest 
component  to  be  considered  Is  the  component  which  lies 
in  the  z'-dlrection,  the  direction  of  the  unit  vector  ar» 
Accordingly,  the  projection  of  (3.1)  onto  ar  , 

v£n(0l‘*i’  “  oos  O'  •  .  0.2) 

will  be  integrated  over  p* . 


When  the  generalized  addition  theorem. 


.  i>n(cos©')  +  ^T[ 


A0cos(^p*  )  +Assin(^p')]  . 
(3.3) 


which  was  given  as  Eq.  ($.1)  of  Ref.  1,  is  substituted  in  the 
right-hand  side  of  (3.2),  and  the  resulting  expression 
integrated  term  by  term  over  the  angle  P*  ,  with  ©,p,9f 
held  constant,  each  term  in  the  summation  over  jx  integrates 
to  zero  (since  Ac  and  Ag  do  not  depend  upon  p * ) ,  and  the 
result  is  the  first  of  the  six  Integral  theorems; 


>'fv 

I  d-a'  “  2Tr^(e,p)  008  0'  Pn(oos»»)  . 


£'<=0 


(**.3b) 
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Prom  (4,1-3)  and  from  the  definitions  of  B®  and  C®  in  (2,4-5) 
it  can  be  seen  that 


■k  S'W 

Vn(n+1)  A  , 

*  -  ar* S  (®t  *,®i)  * 

sin©'  ~r  ~n  l  l 

-  -  Vn(n+1)  ar-C®(01,p1)  . 

(<>.5) 

Differentiation  of  (3.3)  with  respect  to  ©'  now  gives 


Sr sin2®'  ?„(oo3  O') 


d(cos  6* )  n 


>,n. 


V-1  r  0AO  0AS 

-  sin  6*  ^  cos(/up.)  +  sln(^p* ) 


■]}• 

(4.6) 


while  differentiation  with  respeot  to  p*  gives 


/  .n  : 

^ — 1  r 


Zx'Snie  1»%)  =  2^  [^Ac  )  -  >»  Ascos(^p»)]  .  (4.7) 


>*=1 


When  the  explicit  forms  for  A.  and  A  .  obtainable  from  Eq.  (4.8) 

of  Ref.  1,  are  introduced,  then  (4,6)  and  (4.7)  become  addition 

theorems  in  their  own  right,  for  use  with  the  vector  spherical 

harmonics  B?  and  C®  , 

— n  ~n 
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5.  RECURSION  FORMULAS 

The  familiar  recursion*  formula  for  spherical  harmonics^, 


tosej^(s,p) 


(n-jni  +1 ) 
(2n+l) 


+ 


(n-Hml) 

(2n+l) 


(5.1) 


leads  directly  to  a  similar  formula  for  the  veotor  spherical 
harmonics  P°(0,p)s 


oose  p£(e,p) 


(n-|n|+l) 

(2n+l) 


+ 


(n+jml) 

(2n+l) 


(5.2) 


Corresponding  recursion  formulas  can  be  established  for  the 
vector  spherical  harmonics  B°(e,p)  and  C^(e,p),  through 
(5.1)  as  applied  to  (2.^)  and  (2.5) *  but  there  is  some 
cross-coupling  of  the  two  symmetries; 


cose #(e,pj  =  <5±L*?ligr?«L  B“ (e,p)  -  -i“_  tfW) 

~n  F  (2n+l)  (n+1 )  -nvl'  ^  n(n+l)  ~n' 

+  fa.N) /(n-lHn-H)  m 

(2n+l)  n  ~n-1  ^ 


•08  0CfJ(9.;*1  - 

~n  (2n+l)  (n+1)  -n+1' 


i  m 


n(ri+l) 


,  (n+|m|)  </" [n-1 )  (n+1 )  _H 

+  — (S5TI - •  (5^ 
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6,  THEOREM  IV  ; 


Prom  the  equations  In  Section  3  of  Reference  1,  It  can 


be  shown  that 


*  sin©1  cos  (0*  -  0  cos  0 )  *  — •  (cos  6  cos 0*  -  cos  0,  )  ,  (6.1) 
-©  -ir  r  r  sin©  x 


(aA*a1  _)  «  sine*  sln(p'-p  cos  e)  *-*  ~~““-(cosO,)  .  (6.2) 

-P  ~ir  sin©  0p* 


Prom  (2.3),  (5.1).  and  (6.1)  It  is  apparent  that 


5e-ES<ei^i> 


a 

i  e  [ 


x?+i(ei*^i)  "  cose  cose* 

(n+Imj)  m  *| 


(2n+l 


Each  of. the  scalar  spherical  harmonics  on  the  right-hand  side 
of  (6.3)  can  be  replaced  by  an  expansion  of  the  form  (3.3). 
Integration  over  p*  then  gives 


*****  . 

1* 


p'=o 


-  ^cose  x^(e,p)J  £cos  e*  ?n(cose*)^j 


*  i-jj. 
(2n+l ) 


^(e.p)  Pn_1(cos©')\  (6.4) 
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Eq,  (6.4)  oan  be  simplified  with  the  aid  of  (5.1), 
and  with  the  use  of  (4.3a)  in  tho  form  it  takes  with  m=0. 
The  result  is  the  first' half  of  the  fourth  integral  theorem* 


p»  =  2,r 


fae-£‘VW  - 

p»=o 


-<n^>  .m*«- 


P„ foos  ©r 


d(cos  O’ )  n 

(IVa) 


From  (2.3),  (6.2),  and  (4.5),  the  relationship 


2p*£n^°l»^l^ 


1  9  -  0 

sin  0  Qp*  I  1 


-,[cos  e1%ie1,p1)] 

Vntn+l)  r  -  -1 

~i‘  Lcos  el 


sin  © 


(6.5) 


can  be  established,  and  tho  recursion  formulas  (5.1)  and  (5.4) 
can  then  be  used  to  replace  the  bracketed  expressions  by 
linear  combinations  of  scalar  or  veotor  harmonics.  The.  . 
integration  over  p’  can  then  be  carried  out  without  difficulty, 
by  methods  used  in  earlier  theorems,  with  the  result; 

pt  2<v 

Jftp'&WW'  ■  ; ^"ine  •’> 

(IVb) 


•  t 

f: 


p*  =0 


at 

p\ 
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Theorems  (IVa)  and  (IVb)  deal  with  the  transverse 
components  of  P^(©^  f/0-^ ) ,  that  iSj  the  components  that 
are  perpendicular  to  the  unit  vector  ar*  A  vector  combination 
of  Just  these  transverse  components  is  equivalent  to  the 
subtraction  of  the  longitudinal  component  from  the 
original  vectors 


=  2e  5e-*S<9l>V  +  Sp*j£<«L' 


(6.6) 

When  this  vector  combination  of  (IVa)  and  (IVb)  is  made,  it 
is  found  to  have  the  grouping  of  terms  that  appears  in  (2.4). 

The  result  is; 


r  pr  " 

J  -  £r-?r'iS(el>V]  a?' 

p*=o 

=  £p(e,p)  sin2©1  j  -----  P  ( 

Vn(n+1)  ~n  d(oos  ©* )  n 


cos©*).  (IV) 
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?.  THEOREM  V 


The  fifth  integral  theorem  deals  with  the  transverse 
part  of  the  vector  spherical  harmonic  •  *n  analogy 

with  (6.1-2) t  four  scalar  products  will  be  needed: 


tee'Sie’  “ 


^Sp'Sie^ 


sln€^  cos©  sin©*  3©^ 
sin  ©  sin  ©  3©’  * 

sin ©i  cos  6  sin©*  3 pi 
sin  ©  ae*  * 

cos  ©x  a©! 
sin©  ap» 

sin  ©^  cos©*  sin  ©^  cos  ©^  dfii 
sin  ©  sin  ©  ap* 


(7.1a) 

(7.1b) 

(7.1b) 

(7. Id) 


which  can  all  be  derived  from  the  equations  in  Section  3 
of  Reference  1, 

From  (2.4),  (4,3) »  and  (7.1)  it  is  found  that  the 
transverse  portion  of  is* 
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2e  5e*iS<W  +  Sp  Sp'&h’h'1 


[2o 


cos  G  sin  0*  9 


Bin© 


cos  ©i  g  _  * 

+  ~P  sinG  9p*  n^el*^l^ 

«• 


Vn 

(2n+l) 


n(n-H) 

1)  sin  G  ~e 


(n-lnj-fl) 

(n+1 )  '*a#< 


C<,W  * 


*VV 


vmW7t7‘>*W 


(7.2) 


In  the  integration  of  (7.2)  over  p* ,  there  are  five  terms 
to  be  considered.  The  first  two  terms  can  be  transformed 
with  the  aid  of  (4.4),  (4,5),  and  (5.4),  then  integrated 
through  Theorems  II  and  III.  The  remaining  three  terms 
may  be  integrated  directly,  by  Theorem  I  or  its  scalar 
equivalent,  Eq.  (5.2)  of  Eeference  1.  What  is  obtained 
from  these  Integrations  is  the  fifth  integral  theorems 


P’p-rr 

J  -  2r  Sr’jS'W]^’ 

p»=0 


e  2  -7T  E?(e,p)  cos©*  P  (cos©1)  +  ~ —  ■ 6-  -- — — -  P  (cos  6* 

“n  L  n  n(n+l)  d(cos  ©')  a' 

(V) 


-]• 


3 

I 


1 
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8.  THEOREM  VI 


Four  additional  partial  derivatives,  analogous  to 
those  given  in  (^,1)  and  (4,2),  can  be  obtained  from 
Eqs,  (3.5)  of  Reference  Is 


£3l 

l 

sin  0  cos  e’  +  cos  Osin©1  cos(p*-  p  cos  ©) 

ee 

sin  ©l 

-  p  sin2©  sin  G1  sin(p*- p  cos  ©)*j  , 

(8.1a) 

sei 

1 

sin  C  cos  ©  sin  0*  sln(p?-  p  cos  0)  , 

(8.1b) 

Up 

sin  ©1 

a© 

sin©1 

sin2©^ 

J^-  cos©  003©*  sin(p*-p  cos  ©) 

+  sin  ©  sin  ©'  sin(p'  -p  cos  ©)  cos(p*-  p  cos  ©) 


+  p  sin  6  cos  6  sin  ©* 

+  p  sin2-©  cos  ©’  cos(p*-p  cos  ©)”"]  •  (8,1c) 


ep 


sin  © 


r 


sln^e-^  L 


sin©  +  cos  ©sin©*  cos©1  cos(p*-p  cos  6) 


-  sin©  sin^©’  cos^(p*  -  poos  ©) 


]• 


(8, Id) 
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In  terms  of  these  derivatives,  (?•!)  can  be  rewritten  as 


(£©•£ 1©) 


te©*2ip^ 


^~p*£l©^  = 


sln©1  3^1 
sin©  5T  * 

- 1  3©i 

sin  ©  ap  9 

r  W 

=  -  sin  ©,  rr: 


-  p  tan  ©  +  p  tan  © 


0©-,  9©. 

~  +  p  tan©  ~ 

a  a  ' 


0P  * 


(3.2a) 


(8.2b) 


1  /8  , 
(8- 


(8. 2d) 


In  analogy  with  (7.2),  the  transverse  portion  of  the 
vector  spherical  harmonic  oan  be  w1’3-*'*'611  in 


the  form; 


«e  +  Sp  2p*^(e 

“  (Se  iiT®  jy? 


r  e 

-a*  —  - 

-P  L  aa 


imp  tan  8  +  p  tan  © 


<8-3> 
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When  in  (8.3)  is  replaced  by  the  expansion  (3.3) » 
and  the  resulting  expression  integrated  over  p* ,  the  sixth 
of  tho  integral  theorems  is  obtained: 


=  2^r^(©,p)  Pn(cos  e* 


-  SrV£>l^l>]  d?' 


(VI) 
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9.  APPLICATIONS 

These  six  Integral  theorems  are  of  particular  value 

when  a  vector  field  function  Is  to  be  obtained  from  a 

rector  source  function.  An  example  is  the  retarded 

an  electrical  ^ 

Hertz  veotor,  as  obtained  from  /  source-current  distribution. 

The  theorems  show  that  if  a  source  function,  depending  on 

the  angles  has  that  dependence  characterized  by 

A’ 'particular  choice  of  n  and  m,  then  the  integration  over 

the  azimuth  angle  p'  will  Immediately  ensure  that  the  field 

* 

funotion,  in  its  dependence  upon  the  angles  (6, p) ,  will  be 
characterized  by  the  same  choice  of  n  and  m,  provided  that 
the  Green* s  function  relating  field  to  source  does  not  itself 
depend-upon  p'.  Whatever  dependence  the  Green's  function 
may  have  upon  the  angle  G*  (the  angle  between  the  vector  r^  to 
the  source  point  and  the  vector  t  to  the  field  point.)  is  here 
immaterial,  and  cannot  affect  the  'mode  separation*  of  the 
expansion  into  veotor  spherical  harmonics. 

It  should  be  noted,  however,  that  certain  of  the 
six  theorems  introduce  coupling  between  a  P®  source  function 
and  a  E|P  field  function,  or  vice  versa.  That  is,  the 
theorems  mix  the  P®  and  eP  symmetries,  while  keeping  the 
C®  symmetry  separate  from  the  other  two. 

The  use  of  these  theorems  in  electromagnetism  will  be 

h. 

Illustrated  in  a  following  article  . 
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ABSTRACT 

An  equation  giving  the  retarded  Hertz  vector, 

T Tii  ,  t),  in  terms  of  a  source-current  distribution!  JCr^it^)* 
is  derived.  Both  vector  functions  are  written  as  expansions 
in  vector  spherical  harmonics,  with  the  expansion  coefficients 
containing  the  dependence  upon  radial  distance  and  time, 
while  the  angular  dependence  is  kept  within  the  harmonic 
functions.  After  integration  over  two  angles,  expressions 
are  obtained  giving  the  expansion  coefficients  for  the 
Hertz  vector,  which  depend  upon  (r,t),  in  terms  of  the 
corresponding  expansion  coefficients  for  the  current, 
which  depend  upon  (r^,t^).  The  original  four-dimensional 
problem  is  thus  reduced  to  wo  dimensions,  but  with  the 
four-dimensional  causality  requirements  satisfied  at  each 
step  of  the  analysis,  - 
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1.  INTRODUCTION 

t 

Although  the  electromagnetic  radiation  from  a  general 
time-dependent  source- current  distribution  is  usually'  - 
analyzed  in  the  frequency  domain,  for  certain  problems 
tho  time  domain  is  more  appropriate.  An  example  is  the 
electromagnetic  radiation  from  the  electrical  ourrents 
generated  In  the  air  by  an  atmospheric  nucloar  detonation1. 

In  this  example  the  source  currents  are  highly  transient, 
and  retardation  across  tho  source  region  plays  an  important 
rolo.  Furthermore,  part  of  tho  current  is  in  the  form 
of  relativlatio  electrons,  produced  by  Compton  collisions 
between  gamma-rays  from  the  detonation  and  electrons  from 
air  molecules.  The  fields  radiated  by  relativistic 
electrons  show  forward  directivity,  which  is  not  easily 
expressed  in  the  frequency  domain  but  which  enters  readily 
Into  a  time-domain  formalism.  •  la  *.i<.  th:. 

The  analysis  Is  most  straightforward  if  the  medium  is 
assumed  to  be  the  vacuum.  In  tho  example"*  mentioned  above,  there  is 
time-dependent  conductivity  in  the  source  region,  but 
this  conductivity  is  aotually  formed  of  electrons  and  ions 
whose  motion  in  the  local  fields  can  be  represented  as  an 
addition  to  the  primary  source  current.  Similarly, 


dielectrio  polarisation  of  tho  air  moleoures  oan  also  be 
represented  mathematically  as  a  secondary  Increment  to 
the  source-current  distribution.  Thus  if  the  secondary 
currents  are  all  treated  explicitly,  the  vacuum  equations 
will  be  adequate  for  the  calculation  of  the  fields 
generated  by  hll  the  scuroe  currents. 

The  source-current  distribution  will  be  expressed  as 
an  expansion  in  vector  spherical  harmonics,  with  each  coefficient 

depending  in  an  arbitrary  way  upon  the  radial  distance  and 

2 

tho  tine.  With  the  aid  of  theorems  established  earlier  , 
the  angular  dependence  of  the  source  will  be  integrated 
over,  leaving  the  field  function  for  each  vector-spherical- 
harmonic  mode  expressed  as  a  function  of  radial  distanco 
ard  time. 

The  Hertz-vector  formalism  will  be  used,  because  of 
the  relative  simplicity  of  the  integrations.  In  later 

3  it 

articles the  magnetic  field. and  electric  field  will  be 


given  explicitly. 
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2.  HERTZ.  VECTOR  FORMALISM 


The  electromagnetic  fields  generated  by  a  time-dependent 

distribution  of  charges  and  currents  could  bo  expressed 

•* 

in  terms  of  the  vector  and  scalar  potentials,  together 
with  an  auxiliary  condition  which  limited  the  charge  and 
current  densities  to  those  which  satisfied  the  equation 
of  continuity.  It  is  more  convenient,  however,  to  make 
use  of  the  Hertz  vector  formalism,-*  The  need  for  an 
auxiliary  condition  is  avoided  through  the  use  of  the 
free-charge  polarization  vector,  P,  to  represent  the 
source  distribution.  As' will  -be  shown,  the  integrals 
can  then  be  reformulated  so  that  the  final  equation  for 
the  Hertz  vector  expresses  the  source  in  terms  of  the 
current  density,  J,  without  any  explicit  appearance  of 
the  charge  density,  p. 


The  free-charge  polarization  vector  is  defined  through 
the  equations 


P 


V.p  . 


(2.1) 


For  a  transient  source,  therefore, 
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l 


■7« 


dt 


(2.2) 


where  is  a  time  which  precedes  any  current  flow. 

(Any  pre-existing  electrostatic  field  can  be  represented 
through  a  constant  of  integration,  P  ,  added  to  the 
right-hand  side  of  (2.2),)  Because  of  the  form  of  (2,1 
the  equation  of  continuity  of  charge  is  automatically 
satisfied. 

The  Hertz  vector,  rJf1(r,t),  is  defined  by 


1JV’  -  wr  III  I  nix,t-s/c)  dvx , 


(2.3) 


in  the  MKS  units  used  by  Stratton^,  In  (2.3),  dV^  is  the 
volume  element  at  the  source  point,  r^,  and  s  is  the 
distance  from  the  source  point  to  the  field  point,  t; 


s  =  r  -  r,  .  (2.4) 

I  **  I 

Eq.  (2,3)  defines  the  retarded  Hertz  vector,  from  which 
the  retarded  electric  and  magnetic  field  vectors  can  be 
obtained,  through  the  equations 
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(2.5) 

(2.6) 


These  field  vectors  satisfy  Maxwell's  equations. 

An  advanced  Herts  vector  could  also  be  defined, 

through  the  use  of  P(i\  ,t*fs/c)  instead  of  P(r-,  ,t-s/c) 

«  ^  ~-L 

•in  the  integrand  of  (2.3).  Eqs,  (2.5-6)  would  then 
give  advanced  electromagnetic  field  vectors.  These 
would  also  satisfy  Maxwell's  equations.  However,  there 
is  no  evidence  that  these  advanced  fields  play  any 
role  in  macroscopic  electromagnetism.  Accordingly, 
only  the  retarded  Hertz  vector  (2.3)  will  be  considered 
further. in  this  article. 
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3.  EXPANSION  I'  VECTOR  SPHERICAL  HARMONICS 


The  scalar  and  vector  spherical  harmonics  have  been 
defined  in  Section  2  of  Reference  2,  The  vector  functions, 
J,  P,  and  ,  will  be  expanded  in  terms  of  the  vector 

**  +*  ZJ 

spherical  harmonics,  P™,  B°,  and  C?.  The  Hertz  vector 

~n  **i\ 

is  defined  at  the  field  point  (r,t),  and  will  be  written 
as  the  expansion 


co  +n 


U  ■  E  E{TO- 


n=0  m=s-n 


+  Tc.n  <r.t>eg<e.,>^ 


(3.1) 


while  the  free-charge  polarization  vector  is  defined  at  the 
source  point  and  has  the  expansion 


co  "to. 

\ — »  r 


l  -  Y_  £  {?r,n<rl-V  +  Ib.^1^  £<V 


n~0  m=-n 


The  current  density,  J,  has  an  expansion  of  just  the  form  (3.2) 
but  involving  the  expansion  coefficients  » 
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Prom  the  definitions  of  the  vector  spherical  harmonics 
(in  Pief*  2)  it  can  be  seen  that  B®  and  are  both  identically 
zero,  so  that  for  n*0  only  the  terms  involving  Pq  are 
nonvanishing.  The  definitions  also  show  that  the  vector  * 
spherical  harmonics  which  differ  only  in  the  algebraic 
sign  of  m  are  complex  conjugates.  Since  the  vector 
functions  J,  P,  and  7 r  are  real  quantities,  not  complex 
quantities  (in  this  time -domain  analysis),  it  consequently 
follows  that  the  expansion  coefficients  which  differ  only 
In  the  sign  of  m  are  also  complex  conjugates,  as  illustrated  by 
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4,  DEFINITION  OF  O'  AND  p' 

It  will  be  convenient  to  define  a  polar  coordinate 
system  whose  polar  axis  lies  along  the  field-point 
vector  r,  In  terms  of  this  coordinate  frame,  the 
source-point  vector  r^  has  the  polar  coordinates  (r^,6*,p* )* 
and  the  source-point  volume  element  is  given  by 

dVl  ~  TlBin  G’  **1  de*  dp'..  (4.1) 


The  polar  angle  0*  is  the  angle  included  betvreen  the 
vectors  r  and  T-,  ,  so  that 

r^r^  =  r  r^  cos  9*  ,  (4.2) 


and  the  azimuth  angle  p'  will  be  defined  as  shown  in 
Fig.  3  of  Reference  6.  The  trigonometric  transformation 
which  expresses  (e^  ,^5  In  terms  of  (G,p)  and  (©’ ,p* ) 
is  given  as  Eqs.  (3.5)  of  Ref,  6. 


In  this  coordinate  system,  the  Hertz  vector  (2,3)  taker 
the  form 
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5.  INTEGRATION  OVER  p* 

In  the  expression  for  the  Hertz  vector,  Eg..  (4.3), 
the  integration  over  the  azimuth  angle  p*  •will  be 
carried  out  first.  The  geometrical  configuration  is 
shown  in  Pig.  1  of  Reference  2.  The  source-point 
vector,  r,  ,  is  maintained  at  the  constant  length,  r^, 
and  is  swung  about  the  field-point  vector,  r,  with  the 
angle  ©'  between  these  two  vectors  held  constant. 

The  integration  over  p*  is  a  full  circuit,  from 

p*  =  0  to  p*  =  2?r,  and  in  this  integration  the  distance  s, 

defined  in  (2.4),  remains  constant  (since  ©*  is  constant). 

ius  in  the  integrand  of  the  p1 -integration  only  the 

source-point  angles  (^jP^)  will  vary.  By  (3.2)  it  can 
the 

be  seen  that/variation  is  therefore  confined  to  the 
w  tor  spherical  harmonics  themselves,  and  that  the 
expansion  coefficients  will  remain  constant. 

The  theorems  of  Reference  2  can  now  be  brought  into 
the  analysis,  and  utilized  in  the  p' -integration  of  (4,3). 
In  'this  way  all  nine  components  of  the  three  vector 
spherical  harmonics  for  a  given  (n,m)  oan  be  integrated 
over  p* ,  giving  explicit  functions  of  cos  ©*,  multiplying 
vector  spherical  harmonics  of  (©,p). 
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6.  INTEGRATION  OVER  6' 


After  the  integration  over  p*  has  been  done,  the 
integration  over  O'  takes  the  general  form 


Ffr.r^t) 


ii 

|g(cos  0»  ) 


it-s/o)  sin©'  dQ*  , 


(6.1) 


where  f  is  one  of  the  expansion  coefficients  in  (3.2) 
and  g  is  an  explicit  function  obtained  from  one  of 
the  six  integral  theorems  In  Reference  2. 


During  the  Integration  (6.1),  the  distances  r  and  r^ 
are  held  constant,  but  the  distance  s  changes  as  ©'  changes. 
By  the  law  of  cosines,  the  relationship  (2,4)  can  be  written 

as 

s2  *  r2  +  r^2  -  2rr1cosOf,  ->-*(6.2) 

and  its  differentiation,  with  r  and  r^  held  constant,  gives 
s  4s  =  rr^sine*  d©’ .  (6.3) 

Eq,  (6.1)  can  thus  be  replaced  by 

ss£+rj> 

Ftr.r-pt)  «  I  ~  s(E)  .flrjit-s/o)  ds  ,  (6.4) 

8  .|r"ri|  1 
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where  tho  quantity  £  Is  defined  by 


£  =  ‘  cos  6*  = 


2  r  r-, 


(r2  ♦  r* 


s2)  . 


(6.5) 


The  integration  (6.^)  can  be  rewritten  also  as  an 
Integration  over  the  source-point  time  variable,  t1 ,  through 


t,  **  t  -  s/o  , 


-  Ids- 


(6.6) 


(6.?) 


In  this  case  the  quantity  £  is  expressed  as  a  function 


of  t-^  having  the  form  £t,  where 

c'i^p+  ri2  -•8<*-V8]- 

The  integration  then  has  the  appearance 
F^.^.t)  =  J  -6(St)'f(rj 

Vt-o(r+rl)  1 


(6.8) 


l»t^)  dt^  •  (6.9) 


The  function  g(£),  obtained  from  the  integral  theorems 
of  Reference  2j  is  in  each  case  a  polynomial  in  £,  and 
through  (6.5)  it  is  therefore  an  even  polynomial  in  s. 

For  each  such  function,  an  associated  function  G(r,r^,s) 
can  be  defined  by 
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G(r»*i»s)  * 


j 

s-0 


(6.10) 


and  this  associated  function  is  now  an  odd  polynomial  in  s. 


The  integration  over  6*,  as  transformed  by  (6.4)  into 
an  integration  over  s,  will  be  carried  out  as  an  integration 
by  parts,  with  the  aid  of  (6.10).  The  result  is 


Ffrjr^t)  « 


JL_  G(r,r,  ,s)  f(r,,t-s/c) 
r  -  *  * 


s=(r+r-j_ ) 


s=(r+ri) 

{  FFi  G(r,ri,s)  [el f  Wo)]  ds  . 
s=  lr-ri| 


(6.11) 


In  each  case  to  be  considered,  the  function  f (r-^ , t^ ) 
represents  an  expansion  coefficient  vtfiich  has  the 
time  dependence  of  a  component  of  the  free-charge 
polarization  vector,  P,  The  time  derivative  of 
such  a  component  gives  the  corresponding  component  of 
the  source-current  density,  *T,  as  shown  by  (2.1). 

Thus  the  derivative  6f/0s  in  (6.11)  actually  represents 
a  current  component.  An  example  is  given  by 

is  ?r,n^rl*t~s/°J  =  “  c  Jr,n(rl,t“R/c}  • 


(6.12) 
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Similarly,  where  f  appears  in  (6.11)  undifferentiated 
representing  for  example  the  coefficient  P°  n  ,  it  can  be 
written  as  a  tine  integral  of  the  coefficient  * 


IrVrl,Wo)  f 


t, et-s/0 


Jr,n<rX*tl>  dt] 


(6.13) 


tl=to 


where  tQ  is  to  be  chosen  as  a  time  which  precedes  any 
current  flow.  (As  mentioned  in  connection  with  Eq.  (2.2), 
any  pre-existing  static  field,  formed  by  an  earlier 
charge  displacement  which  is  not  included  in  the 
transient  current  that  is  being  analyzed,  can  be 
represented  by  a  constant  of  integration,  added  to 
the  right-hand  side  of  (6.13).) 
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7.  MATRIX  REPRESENTATION 


After  the  integration  over  p' ,  and  the  transformation 
from  (6.1)  to  (6,4),  but  before  the  integration  over  s, 
the  Hertz-vector  expansion  coefficients  can  be  represented 
oompactly  through  the  matrix  equation: 


IT- 


*iV° 

h  A1 

ri=° 


s«(r+r1) 

/r)d XjJ  dE  • 


Hr-ri| 


(7.1) 


In  (7.1)  each  of  the  Greek  indices,  tf  and  A,  runs  through 
the  three  values;  r,  B,  C,  A  repeated  Greek  index  indicates 
summation  over  these  values.  Thus  (7,1)  represents  three 
equations,  each  of  which  may  involve  three  expansion 
•coefficients  for  the  free-charge  polarization  vector  P. 

However,  some  of  the  **.  elements  of  the  matrix  4“x(S) 
are  zero,  so  that  the  equations  are  in  fact  relatively  simple. 
The  matrix  elements,  determined  through  the  use  of  the 
theorems  in  Reference  2 »  are; 


s£“£(t)  »  5  p„(S)  . 

(7.2a) 

=  s  pnuv  +  n(B^)dtys)-. 

(7.2b) 

\ 

«o“o<£)  =  pn<£>  > 

(7.2o) 

•f 

k. 
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=  =  (^LL  d  P  (5)  , 

B*  Vn(n+1)  <*£  n 

(7. 2d) 

«i>> 

=  e£lU)  =  o  , 

(7.2e) 

g^>(5) 

=  ®C,B^  =  G  • 

(7.2f) 

When  the  integrals  having  the  fora  (6.4)  are 
integrated  by  parts,  as  illustrated  in  (6.11),  then 
the  Hertz-vector  expansion  coefficients  (7.1)  take  the 


forms 


Tl%&  *i»t-  (r+r^/c 

"5? ~  J  (rl/r’arl  Gi?i<r>rl*r+Jl>  7  OW  dtl 

rl“°  *1**0 


ri .?  hjf-  tr-ri)/o 

Fe~  J  (r±/'r'>Ar:  Gff,l(r>ri>r-r i>  J  J\,n^Ti’*i^  dfci 

°rl=0  tl=to 


i  *1^°°  ti^t-Cri-rJ/c 

2^  J<V'r>drl  0an,lir’Tl,Tl~r)  J  JZn<rl’h>  dil 

rl=r  *1^0 


riy  t^Mr-r^/c 

2  eo  J  (r]/r)dri  7  Gi°x(r,rl’st"ctl>  dfci 

T  «s  i  V  .  J  .  rn 


rls° 


rI=c0 


tjL=t-  (r+r-jJ/c 


1  y  ti=t-(ri-r)/c 

2«o  J  <rl/r,drl  J  G^’(r-rl-ot-ctl>  dtl  •  <7* 

^^t-  (r-rrx )/ c 


rl~r 
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and  can  therefore  be  obtained  explicitly  with  the 
aid  of  (7.2). 
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8.  SUMMARY 

Equation  (7.3) *  together  with  (3.1)  and  (3.2),  gives 
the  Hertz  vector  ||(r,t)  in  terns  of  the  source-current 

0<J  ~ 

density  function  *  11118  use  of  the  exPanslon 

in  terms  of  vector  spherical  harmonics,  for  both  the 
Hertz  vector  and  the  source-current  density,  has  provided 
a  mode  separation, ■ in' which  a  particular  source  mode, 
characterized  by  (n,m),  leads  to  a  field  mode  which  is 
also  characterized  by  (n,m).  A  source  current  with  the 
symmetry  of  the  vector  spherical  harmonic  leads  to 

a  field  mode  with  this  same  symmetry,  but  there  is 
cross-coupling  between  source  and  field  inodes  having 
the  symmetries  of  the  vector  spherical  harmonics  and  Er^. 

The  retarded  Hertz  vector,  r * t ,  is  expressed 
as  an  integral  over  the  source  current,  £(£1*^1) * 
within  the  region  of  space-time  which  is  consistent  with 
causality  requirements.  Since  the  original  form  of  the 
Hertz  vector,  given  in  Eq,  (4.3),  is  in  accord  with 
the  requirements  of  causality,  and  since  there  is  no 
contamination  by  ’advanced*  fields,  either  in  (4,3)  or 
In  any  later  stages  of  the  analysis,  it  can  be  concluded  that 
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the  Hertz  vector  in  the  form  (7.3) »  though  expressed  in 
the  (r,t) -plane  where  causality  requirements  are  not 
very  transparent,  will  nevertheless  remain  fully 
consistent  with  the  physical  requirement  that  the 
electromagnetic  effect  of  a  moving  charged  particle  should 
travel .at  the  velocity  of  light  if  the  medium  is  the. vacuum 

The  explicit  calculation  of  the  electric  and -magnetic 

fields, -’from' the  Hertz  vector  given  here  in  (7*3)  and  (3*1) 

u 

will  be  carried  out  in  two  following  articles.  * 
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ABSTRACT 

The  preceding  article  gave  tne  retarded  Hertz  vector 
in  terms  of  a  general  transient  source-current  distribution. 
Tb ;  •Hart'”  vector  is  here  differentiated  to  give  the  vector 
potential  and  the  magnetic  field  vector.  All  of  these 
vector  quantities  have  been  expanded  in  vector  spherical 
harmonics,  and  in  each  case  it  is  the  expansion  coefficients, 
which  are  functions  of  the  radial  distance  r  and  the 
time  t,  whion  are  expressed  as  two-dimensional  integrals, 
in  the  (r,t'-plane,  over  the  source-current  expansion 
coefficients  ha’  .g  the  corresponding  values  of  n  and  m, 
where  these  are  the  mode  parameters  characterizing  the 
vector  spherical  harmonics,  P™(6,p>,  |P(e,p),  and  C^(6,p). 
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1.  INTRODUCTION 

In  the  previous  article1  the  retarded  Hertz  vector  was 
obtained,  for  a  fully  general  transient  source-current 
distribution.  The  source  current ,  J(rn ,t, ),  was  expressed  as  an 
expansion  in  vector  spherical  harmonics,  with  the  dependence 
upon  the  angles  (0^,)^)  contained  within  these  harmonics, 
while  the  dependence  upon  the  radial  distance  r^  and 
the  source  time  t^  was  contained  within  the  expansion 
coefficients.  Similarly,  the  Hertz  vector, r[fT(r  ,t) ,  was  expressed 
as  an  expansion  in  vector  spherical  harmonics  which  were 
functions  of  the  angles  (0,p) ,  with  expansion  coefficients 
which  were  functions  of  the  radial  distance  r  and  the 
time  t. 

It  was  found  that,  for  a  given  (m,n),  labeling  the  vector 

spherical  harmonics,  the  expansion  coefficients  for  the 

Hertz  vector  could  be  expressed  in  terms  of  the  expansion  . 

coefficients  for  the  source  current.  There  was  no  coupling 

between  the  expansion  coefficients  for  different  choices  of  n 

or  for  different  choices  of  m,  There  was,  however,  coupling 

between  two  of  the  three  vector  harmonics,  those  denoted 

by  Pra  and  JeP  ,  but  no  cross-coupling  between  either  of  these 
~n  ~n 

two  and  the  third  vector  harmonic,  denoted  by  . 
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In  the  present  article,  this  Hertz  vector  will 
be  substituted  into  the  equation 


eo  ix  is IT 

rsj 


(l.i 


which  gives  the  magnetic  field  H  in  terms  of  the 

Hertz  vector  rfp,  in  MKS  units  for  which 
Aj 

*  8.834  •  10~12  farad/meter  .  (1.5 

An  alternative  formulation  will  also  be  given,  in  terms 
of  the  vector  potential,  A,  defined  here  by  the  equation 


i  ± 

o2  at 


TT  * 


For  the  quasi-vacuum  conditions  that  have  been  postulated 
the  magnetic  field  g  is  then  given  by 
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2,  MAGNETIC  FIELD  EXPANSION 


If  the  Hertz  vector  and  the  magnetic  field  vector 
are  both  expressed  in  polar-coordinate  components!  then 
the  vector  equation  (1.1)  separates  into  the  component 


equations : 


*  £°  /  .  £.  A TT 

r  sin  6  1  sp  et  0 


+  (2,la) 

-  (z-ib) 


However,  it  will  be  more  convenient,  and  will  lead 
to  field  expressions  which  are  more  compact,  if  the 
magnetic  field  vector  is  first  expanded  in  terms  of 
the  vector  spherical  harmonics  defined  in  Section  2 
of  Reference  2  •  This  expansion  can  be  written  in 
the  forms 


co  +n 


E(r,t)  =  £  £{HrVr,t)  + 

n=l  m=-n L  „  *\ 

0^(e,^)  t.  .  (2.2) 


+  Hc“n(r,t) 
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This  summation  begins  with  n*l  because  the  spherically 
symmetric  current  component,  with  n a  0,  does  not  give 
rise  to  any  component  of  magnet io  field,  but  only  to 
a  spherically  symmetric,  radially-directed  electric  field. 


The  expansion  (2.2)  parallels  the  similar  expansions 
for  the  Hertz  vector  and  for  the  current-density  vector, 
as  described  in  Section  3  of  Reference  1.  As  in  the 
earlier  expansions,  the  reality  of  the  field  vector 
requires  that  expansion  coefficients  which  dlffei  only 
In  the  sign  of  m  should  be  complex  conjugates  of  each 
other,  as  in  the  example; 


Or<t) 


(2.3) 


/ 
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3.  time  derivative  of  hertz  vector 


When  the  Hertz  vector,  in  the  matrix  form  in 
Equation  (7.1)  of  Reference  1,  is  differentiated  with 
respect  to  the  observation  time  t,  the  result  is 


r^«co  s=  (r-Kr-^ ) 

HTcr!n tr-fc’  =  air-  J <Vr >  dli  J  ( ? >  JjT. n (rl • t- 

°rl'’C  Ec  |r-rl| 


_0 

et 


•s/c)  ds 

(3.1) 


The  notation  here  is  the  same  as  the  notation  in  Reference  1, 

with  Greek  subscripts  running  through  the  three  values  r,  B,  C, 

associated  with  the  three  vector  spherical  harmonics,  P?,  #,  0°. 

~n  ~n  ~n 

As  in  (2.2) i  the  dependence  upon  time  and  radial  distance 
lies  with  the  expansion  coefficients,  so  that  when  the 
time  derivative  is  taken  it  is  only  necessary  to  differentiate 
these  expansion  coefficients,  as  was  done  here  in  (3.1). 


When  the  'integration'  over  ds  in  (3.1)  is  replaced 
by  an  integration  over  dt^,  with  the  use  of  (6,4-9)  of 
Reference  1,  then  (3.1)  takes  the  form 


ITm"-1’  =  FF  <">  (Vr> 


J”,n(rl-V  dti  dr: 


(3.2) 
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where  the  abbreviation  US)  has  the  following  equivalent 


meanings; 

US)  - 

rl?c0  tlp-lr-rli/0 

J  ] 

r^*0  t^=t- (r+r1 )/c 

(3.3a) 

US)  * 

t,*t  r,  *r+c(t-tn ) 

V  T 

J  J 

(3.3b) 

• 

tj=-co  rj_»|r-o(t-tj)| 

(//)  - 

rl?.r  q-t-tr-r^/o 

J  J 

i*l*0  ^^t-Cr+r^)/© 

rl?" 

/  /  ,  (3.30 

r^*r  t£=t-(r+r^)/ c 

• 

t1*t-r/c  ^1*r+c(t-t1) 

t,  =  t  r-,*r-fc  (t-tn ) 

1  T  A 

(//}  - 

J  J 

*  J  J  '  (3' 

t^  -oo  r^-r+c  (t-t1) 

t^t-r/c  r^r-cft-t^) 

which  define  the  realm  of  integration  in  the  (r^ , ^ ) -plane . 
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4.  VECTOR  IDENTITIES 


3 

While  Morse  and  Feshbachy  have  given  a  number  of 
mathematical  identities  Involving  vector  spherical 
harmonics,  these  are  not  in  the  form  that  is  needed 
here.  For  the  present  application,  there  are  three 


curl 

identities. 

VxP®(0,pl  = 
~  ~n  '  ^ 

^Fogfe.p,  . 

(4.1) 

V* 

-  r  5n(e*^’  * 

(4.2) 

vxog(e.^  - 

r  +  ntr  11  « 

(4.3) 

which  can  readily  be  generalized  to  the  situation  in 
which  each  vector  spherical  harmonic  is  multiplied  by 
a  scalar  function  of  the  form  Q(r,t) .  The  generalized 
identities  are 


vx[«  0  '  *  FF  «  -n  *  . 

vx[ft  Eg]  »  -  ^f(rQ)  Cg  ,  (4.5) 

SX[Q  Cg]  g.  (4.6, 


Pa  se  103 


5.  MAGNETIC  FIELD 


In  the  interest  of  a  compact  notation,  t4ie  three 
vector  spherical  harmonics  will  be  represented 
collectively  by  vj^n(©,p) ,  where  the  Greek  index  tf 
takes  on  the  three  values  r,  B,  C,  as  in  (3.1). 

Ihe  explicit  definitions  ares 


«.  m  __  pin 

~r,n  ~  ~n  * 

fn  • 


V®  •= 
~B,n 


vm 

-C,n 


c®  • 

*»n 


In  terms  of  this  new  notation,  and  with  the  summation 
convention  for  repeated  Greek  Indices,  (2.2)  can  be 
written  compactly  as; 

to  +n 

Hlr.t)  -  « 

n=l  m=~n 

with  similar  expressions  for  the  Hertz  vector  and  its 
time  derivative. 


(5.1a) 

(5.1b) 

(5.1o) 


(5.2) 
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In  particular,  the  time  derivative  of  the  Hertz  vector 
can  be  given  in  the  form 


-  ir  I  3?.n  <//»  dt^.  (5.3 


Substitution  into  (1,1),  with  the  use  of  gives; 


H£n(r,t)  =  /n(n+l)  <//)  Pn(?t)  rx  dtx  di^  ,  (5.4) 


Hr  (r,t)  =  ~ 

B,n'  2r  Or 


|^(//)  pnUt)  ^c>n(ri,t1)  d^dr^  i  (5.5) 


=  £  ^Tn-af  (//>  Pn(?tl  Jr*n(rl’V  dtx 


-  ?r  Pn^t^  ^  [rl  JaU*rl#tl^]  dtl 


dr,  .  (5.6) 


In  Equation  (5,6) »  the  second  expression  on  the  right-hand 
side  was  obtained  through  an  integration  by  parts,  in  which 
the  form  (3,3d)  was  used  for  the  integral  operator. 
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6.  DISCUSSION 

Equations  (5.^-6)  give  the  magnetic  field  expansion 
coefficients,  which  are  to  be  Inserted  into  (jj.2).  I*1 

can  be  seen  that  there  is  an  element  of  symmetry  in  the 
dependence  of  the  three  magnetic-field  components  upon 
the  three  current  components ,  but  that  the  symmetry 
Is  not  as  conveniently  expressed  in  matrix  form  as  was 
the  case  for  the  Hertz  vector  itself,  as  given  in 
Eq.  (7.1)  of  Reference  1. 

While  the  forms  (5*5)  and  (5.6)  were  chosen  here 
because  of  their  compactness,  there  are  other  forms 
for  these  magnetic  field  components  which  avoid  the 
use  of  the  operators  3/3r  and  3/3r^.  These  other 
forms  can  be  obtained  from  (£.5)  and  (5.6)  through 
integration  by  parts  and  through  the  carrying-out  of 
indicated  differentiations. 
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ABSTRACT 

The  retarded  Hertz  vector,  obtained  earlier,  is 
differentiated  to  give  the  scalar  potential  and  the 

electric  field,  vector,  associated  vrith  a  general  transient 

* 

source-current  distribution.  The  scalar  and  vector 
quantities  are  expanded  in  terms  of  scalar  and  vector 
spherical  harmonics.  For  each  mode,  characterized  by 
specific  values  for  n  and  m  in  the  expansions, 
the  scalar  potential  and  electric  field  vector  are 
expressed  as  two-dimensional  integrals,  in  the  (r,t) -plane, 
over  the  causally-accesslble  portion  of  the  source-current 
distribution. 
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1.  INTRODUCTION 


In  the  preceding  article  ,  expressions  were  obtained 

for  the  magnetic  field  components  associated  with  a  general 

transient  current  distribution.  Similar  expressions  will 

now  bo  derived  for  the  electric  field  components,  and 

as  before  these  will  be  derived  from  the  previously  derived 
ths  components  of  2 

expressions  for/the  retarded  Hertz  vector  . 

In  terms  of  the  retarded  Hertz  vector  ^  ,  the  electric 

r * 

field  vector  is  given  by  the  equation 

•  E  -  ss-T  -  h  ^TT  .  u.D 

rJ  C  Ot  r* 

In  terms  of  the  more  familiar  scalar  and  vector  potentials, 
the  eieotrlo  vector  is  . 


(1.2) 


The  Identification  of  the  vector  potential, 

“  “  ^2  StU  ’  .  (1‘'3) 

has  already  been  made,  in  (1.3)  of  Reference  1.  A  corresponding 
identification  of  the  scalar  potential,  in  terms  of  the 
Hertz  vector,  is; 


-  s-TT 

A 9 


«s 


9 


a.'o 
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2.  SCALAR  IDENTITIES 

In  analogy  with  the  vector  identities  given  in 
Section  4  of  Reference  1,  there  are  certain  scalar 
Identities  which  arise  from  the  action  of  the  divergence 
operator  upon  the  vector  spherical  harmonics,  and  upon 
products  of  a  scalar  function,  Q(r,t),  aud  the  vector 
spherical  harmonics. 

The  three  basic  divergence  identities  are 

v.p^o.ja)  =  |  x£(e,>si  ,  (2.1) 

=  -  — ^±!I  xJ|Ce,0)  .  (2.2) 

v-c®(e,p)  *  o  .  (2.3) 

The  scalar  spherical  harmonics,  and  the  vector 

spherical  harmonics,  P®,  E^,  and  C^,  have  been 

defined  in  Section  2  of  Reference  3. 

The  more  general  divergence  identities,  in  which 
each  vector  spherical  harmonic  is  multiplied  by  the 
scalar  function  Q(r,t),  ares 
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Xn  ’  (2-'-l) 

s*[«S]  -  <*r's> 

S’[«£n]  =  °*  <2-«’ 

As  will  be  seen,  it  Is  the  vanishing  of  the  divergence 
in  (2.6)  that  can  be  associated  with  the  fact  that 
the  source-current  components  involving  the  vector 
spherical  harmonics  C®  do  not  lead  to  any  charge 
accumulation,  hence  these  currents’ do' not  give  any 
electrostatic  contribution  to  E(r,t).  These  currents 
are  circulatory  in  character,  and  make  only  inductive 
and  radiative  field  contributions. 
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3.  DIVERGENCE  0?  HERTZ  VECTOR 

*  / 

In  the  evaluation  of  the  divergence  of  the  Hertz 

vector,  one  of  the  functions  defined  in  Section  7  of 
Reference  2  appears  repeatedly.  This  is  the  function 

s*--s 

Gc^c(r*rl,s>  “  |  Pn^  ds  »  (3.1) 

s=0 

where  . 

I  -  ~ -  (r?  +  r,?  -  s?)  .  (3.21 

2rrl  1 

Special  values  of  this  function  are 

G^”J,(r  .r^r+r^)  =  (2n+l)"1  (t*4"1^"*1  +  r“nr^+1),  (3.3a) 

Gc”c(r»r1»r"r1*  “  (Zn+l)”1  (rn+1r^n  -  r“nr3n+13 ,  (3.3b> 

GC^C^r,rl*rl"r^  *  (Sn+l^"1  ( -  rn+1  r~n  +  r~n  r^4*1 ) .  (3.3c) 

Equations  (3.3)  are  established,  in  the  Appendix. 


I  <5 
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The  divergence  of  the  Hertz  vector  does  not  involve 
the  source  currents  with  the  jiV  symmetry,  as  noted 

w  |T1 

earlier  in  connection  with  Eq,  (2.6K  However,  both  the 
J^n  and  the  Jg  current  components  make  contributions 
to  this  divergence.  The  explicit  form  for  this  divergence  is; 

Tlr^  ti«t-  (r+r^)/ c 

~  I  J;1  GC?C(r,rl,r+rl  J  J  Jr,n*ri’V  dtidrl 

l  rl*°  tj«to 

r!fr  r  ^(r-r^/c 

+  °o"o(r>rl*r"rl’J  J  Jr,n(rl>tl'dtldrl 

rlr'°  *1**0 

t^t- (r-^-r  1/ c 

+  J  1 1  t  GC?C^r,rl*rl"r)l  J  Jr1,n^rl,tl^  dt]  drl 

.  rl"r  L 

(//)  Gc”c^r,rl*ot"'ctl5J  mn(rl,tl)  dtl  drl 


pl*®  t1=t-(r+r1)/c 

-  Sn ( n  +1 )  j  GQjQ(r»ri>r+r^^  I  dB,n^rl*^l^  d^l  drl 

r,  -C  tn*=t 

j.  l  o 
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r,  =x 


np*  t^t-tr-^1/0 

J  4”o<r>ri>r-ri)  J  JB,n(ri*V  dti  dri 


rl=° 


r,  *00  vlpl'"'Ar*  u  w 

4i'(n+iy  J  Gc^C(r,rl,rl“r)  J  JB,n(rl,tl)  dtl 

rx«r  t1«t0 


Vfco 


dr. 


-  /n(n+l)  (//)  Gc^c(r>ri»°^“0^i^  ^B,n^rl*^l^  dtl  ^*1  /  2£ 


The  integral  symbol  (//)  has  been  defined  in  Bq.  (3.3) 
of  Reference  1, 


t 


-  ^(6. P). 
(3  M 


fHiwSwwT 
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4.  ELECTRIC  FIELD 

Tj.  r  electric  field  vector  is  obtained  from  Eq.  (1.1). 
This  includes  an  operation  in  which  the  gradient  of  the 
divergence  of  the  Hertz  vector  is  taken,  and  for  this 
operation  it  is  convenient  to  make  use  of  the  identity 

x£(e,p)]  =  Fj^  p”<0,;5) 

+  /nCHTlT  Q(r,t)  j£(e,<?)  .  (4.1) 


The  second  time  derivative  of  the  Hertz  vector  also  appears, 
and  this  is  conveniently  found  from  the  time  derivative 
of  the  expression  given  in  Eq,  (3,2)  of  Reference  1. 

After  an  integration  by  parts,  this  second  time  derivative 
Is  found  to  be  expressible  in  the  compact  matrix  form; 


u/i  (V"  [ifrOri'*i>]“i"i 


(4.?) 


As  before,  a  repeated  Greek  index  indicates  summation  over 

/_! 

the  values  r,  B,  C,  and  the  functions  are  as 

defined  in  (7.2)  of  Reference  2, 
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When  the  operations  Indicated  in  Eq.  (1.1)  are 
carried  out,  the  resulting  expressions,  for  the  expansion 
coefficients  for  the  electric  field  vector,  are; 


Er">n(r,tl  = 


t,=t 

"  "^O  <f  J*»n*r,tl*  dtl 

Hmto 


r, ^r  t,=t-(r-r, )/c 


.  'tv  (n+1 ) 

+  T*2n+1)  66 


^n+1  dti  dri 


r,  =o  t,  =t 
1  1  o 


n  (n-KH 
(2n+l)  6, 


rl?“ 


r"-1  ^-n  dtj  dr. 


rl“r  V*  o 


n  (n+1)  /  ,  ,_-n-2  _  n+1  ,  „n-l  _  -n.  Tm  ^ 

2(2n+i")  (r  rl  +  r  rl  }  Jr , n 'P1  * tl '  dtl  drl 


*  TgS  (//)  Gc"o(r>ri'st-Otl1  dtl  drl 


rx*r  t1=t-(r-r1)/c 

*  nfcflf  J  T  '.VvV«h", 

f  r=0  t-,  =t 

1  1  o 


Page  118 


r.=co 
1/ 


2€rtcr 


(-l)n  ^  «rx 


rl=° 


rlp.r 

— i —  r,  j|'n(r1.t-[r-rl]/o)  drX 

2€0or  J  1  3, n  1 


V° 


V® 


2^7  J  rl  JBin(rllt'trl'rJ/c)  drl 

rl“r 


rl7?  t1?t-(r-ri)/0 


n  (n-M) 


(2n*M)  £ 


o  J 

i 


-n-2rn+l  dtx  drx 


r-, -0 


r1=co  t1=;t-(r1-r)/o 

*  fr-*L.  f  (  r*1"1  r{n  V  dtl  drl 

(2n+l)  %  J  J 

r,  =r  tn  =t 


n 


1  o 


,  ntn*),._  {//,  (r-n-2  r  n+l+ rn-l  ri-n,  j-^.V  «x  drx 

-2(2n+l)  £Q 

*  2^  (//) I”1  ^10°>’ri’0t'0tlOJB”n(ri'tl,atl 


(4. It) 
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m 


C»n 


(r,t)  =  - 


(//)  <Vr)  W[a!^  dtlfirl 


(^.5) 
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5,  DISCUSSION 

Equations  (4.3-5)  give  the  components  of  the 

electric  field,  in  terns  of  the  components  of  the  source 

current.  These  equations  are  hot  given  in  the  general  matrix 

for  the  present 

form,  since  their  relative  complexity  precludes/a  simple 
matrix  representation.  However,  there  are  elements  of 
symmetry  involved,  and  a  compact  matrix  representation 
may  eventually  be  devised. 

The  electric  field  components  are  here  given  as 
explicit  integrations  over  the  source-current  components. 

For  the  inductive  and  radiative  contributions,  the  integration 
is  limited  to  the  realm  specified  by  the  notation  US). 

This  realm  limits  the  source  currents  which  are  ‘visible* 
at  the  observation  point.  However,  in  addition  to  these 
immediately-sensed  contributions,  there  are  electric  field 
contributions  which  are  electrostatic  and  are  associated 
with  earlier  flow  of  current,  establishing  a  distribution 
of  electric  charge  which  generates  an  electrostatic  field. 

The  contributions  to  the  electric  field  which  are  associated 
with  this  Electrostatic  dipble-uoment  distribution  are 
expressed  in  the  integr  Is  which  have  as  their  lower  limit 
ti -  tQ,  where  tQ  is  a  time  preceding  any  of  the  transient 
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ourront  flow  which  contributes  to  the  fields  that  are 
being  calculated. 

The  Important  aspect  of  Equations  (4,3-5)  Is  their 
expression  of  the  causal  relationship  between  the  source 
current,  J,  and  the  electric  field,  E,  which  Is  generated 
by  this  source  current.  When  these  equations  are  used  as 

the  basis  for  a  numerical  solution  of  Maxwell's  equations, 

4 

the  requirements  of  causality  will  be  met. 
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APPENDIX 

For  the  derivation  of  Equations  (3.3),  it  will 
first  be  noted  that  the  Legendre  polynomial  Pn(?)» 
which  can  be  written  in  the  form 


,n 


Pn(0 


2n  nl 


^■(?2-l)n, 


at 


n 


(A.l) 


is  an  odd  polynomial  in  t 
polynomial  if  n  is  even, 

rn(-?>  =  (-i'n  p„(0  . 


if  n  is  odd,  an  even 
so  that 


(A, 2) 


When  the  argument  t  is  written  out  as 


1  "  /  ~2  .  _  2  _2% 
—  (r  +  r,  -  s  )  , 


2-rr, 


T41 


(A.3) 


then  the  Legendre  polynomial  P  (£1  has  the  fora 


(?  rij) 


n  n 


(-2,  r,?,  s? )  , 


(A.  4) 


where  the  function  f  is  a  finite  polynomial  in  its 
three  arguments. 


Pago  113 


It  follows  that  the  integral 

Pn(r»r1*R)  * 

is  an  odd  polynomial  in  R,  so  that 
Fn(r,r1,-H)  -  -  Fn(r ,r^,R)  . 

It  can  also  bo  shown,  from  the  form  of  (A,4),  that 


s=R 

1 


pn(?)  as 


(A.5) 


(A.  6) 


Pn(-r,r1,R)  - 

(-1)"  Fn(r,r1,R)  , 

•  (A.?a) 

Fn(r>"rl>p0  85 

(-1)"  Fjjlr.rj.Rl  , 

(A, 7b) 

Fn(-r»-rl»R>  « 

Fn<r,iVR}  . 

(A.7c) 

The  three 

integrals  to  be  evaluated  are 

Pn*r,rl,r+rl^ 

s=r+rx 

-  I  Pn(p  ds  , 

s~0 

• 

(A. 8a) 

•  •  \ 

» 

Pn(r,ri,r-ril 

Srr-Tl 

«  I  Fn(P  as 

8=0 

(A. 8b) 

Vr'rl»rlTp’ 

S^l"r 

-  j  V$)  ds  • 

8=0 

(A.8o) 
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All  three  are  defined  here  with  no  restrictions  on  the 
relative  magnitudes  of  r  and  r^..  I^rom  (A, 6)  It  can 
be  seen  that 


Pn(r,rT  ,r-rn)  =  -  ^(r.r-,  ,rn -r)  , 


n'  *1*1 


and  from  the  symmetrical  way  that  r  ancl 
£,  In  (A. 3),  It  Is  also  clear  that 

Fn(rl’r*rl-r’  M  '  f,n{^■rl>^-rl,  * 
Fn(rl'r>rl+*'1  "  Fn(r,rl,r"irl^  • 

Equation  (A. 3'  can  be  solved  for  s, 

?  9  1/2 

s  *  (r£  -  2  r  rx  g  +  ly  1 

When  r  and.  r^  are  both  held  constant, 
ds  and  dg  will  be  related  by 

s  ds  »  -  r  dg  , 

If  is  smaller  than  r,  the  expansion 


p  p  ^ 

I  *=  Po<?’  +  T  Pl<*'  +  -~Z  P2(5’ 


oonverges,  and  it  is  possible  to  write 


(A. 9&) 

i 

r^  enter  into 

(A. 9b) 
(A.9o) 

giving 

.  .(A.10) 

the  differentials 

(A. 11) 

+  ...  (A.12> 
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s^r+r^ 

?n(r,rl,r+1V  “  Fn(r,r1,r-r1)  =  J*  Pn(£)  ds 


s=r-r- 


'p1 

-  s'lrrlPn 

{=+1 


(?)  as  =  r 


f1 

1  pn 


n 


(?)  |  a?  =  rx  ~  .  (A 

s  x  r11  (2n*l) 


In  (A, 13)  the  series  (A. 12)  has  been  substituted,  and  the 
orthogonality  properties  of  the  Legendre  polynomials  have 
been  utilised. 


If,  on  the  other  hand,  is  greater  than  r,  then  the 
convergent  series  expansion  is 

r  ^ 

~  »  P0(?1  +  ~  Px(?)  +  -T  P2(?)  +  ...  ,  (A.lb) 

1  ri 


and  the  resulting  equation,  analogous  to  (A.13)»  is 


n 


Fn(r,r1,r-:-r1)  -  F„(r ,^,^-r)  =  r  fg 


rx"  (2n-KL) 


(A.15) 


From  (A. 4) *and  (A, 5)  it  Is  apparent  that  the 
multiplication  of  Fn(r,r^,R)  by  (2rr^)n  will  cancel  the 
factors  in  the  denominator  and  give  a  finite  polynomial. 
In  particular,  if  Equation  (A. 13)  is  multiplied  on  both 
sides  by  (2rr^)n,  it  will  become  a  relationship  between 
the  sum  of  two  finite  polynomials  (on  the  left)  and  a 
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monomial  (on  the  right).  Evidently  all  of  the  terms  in 
the  two  polynomials  must  subtract  out,  except  for  the 
terms  which  add  to  give  tb'*'  monomial  on  the  right.  Moreover, 
this  algebraic  relationship  involving  finite  polynomials, 
while  established  for  r^  less  than  r,  can  be  continued 
analytically  to  include  the  regions  where  r-j  is  greater  than  r. 
Similarly,  Equation  (A. 15) 9  while  established  for  r^  greater 
than  r,  is  a  simple  algebraic  relationship  which  can  be 
continued  analytically  to  the  realm  where  r^  is  lees  than  r. 

When  (A. 9a)  is  Inserted  into  (A.15>*  the  result  is 


Pn*r,rl,r+rl)  +  Pn*r,r1,r~rl)  “ 


r11 

«v*  _ _ _ 

**in 


_ 2j_ 

(2n+H 


(A.16) 


Now,  from  (A. 13),  (A. 16),  and  (A, 9a),  the  desired  results 
are  obtai'  ed; 


Pn(r,rl  »r+ri^  = 

(2H+1)"1 

is 

(rn+1  r  j‘a  + 

r-nr1n+1),- 

(A.l?a) 

Pn(r»rl»r“rl^  = 

(2n-M)-1 

(rn+1  rxn  - 

r-nrxn+1)  , 

(A.  l?b) 

Pn(r.ri,ri-r)  = 

(an-M)”1 

(-  rn+1ri-n 

+  r"n  . 

(A.  l?o) 

These  results  are  equivalent  to  Equations  (3*3),  since  It 
is  apparent  from  the  definition  of  F  in  (A. 5)  that 

Fn(r,r1,R)  =  G^Cr.r^R)  .  (A.18) 
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ABSTRACT 

Maxwell* s  equations  In  differential  form  do  not 

distinguish  between  advanced  and  retarded  solutions. 

Unless  speoial  precautions  are  taken,  a -point-by-point 

numerical  integration,  using  a  finite-difference  version 

of  Maxwell's  equations,  will  lead  to  a  mixture  of 

advanced  and  retarded  fields,  which  is  inadmissible 
•  • 

on  physical  grounds.  The  causality  requirement  can  be 
met  if  Maxwell's  equations  are  written  in  integral-equation 
form,  with  retardation  incorporated  in  all  the  integrands, 
.A  solution  using  numerical  integration  then  will  be 
physically  acceptable.  The  distinction  between  causal 
and  non-causal  solutions  is  illustrated  by  an  example 
in  which  the  problem  symmetry  permits  separation  of 
variables  and  the  reduction  of  "the  four-dimensional 
space-time  problem  to  a  two-dimensional  problem  involving 
only  radial  distance  r  and  time  t. 
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INTRODUCTION 

Maxwell’s  equations,  as  they  stand,  admit  of  two 
kinds  of  exact  solutions  for  the  electromagnetic  fields 
generated  by  a  specified  tine-dependent  distribution  of 
electric  charge  and  current.  These  are  the  retarded 
solution  and  the  advanced  solution.  An  example  of 
a  retarded  solution  is  the  Lienard-Wiechert  solution  for 
the  fields  generated  by  a  moving  charged  particle 
The  fields  at  an  observation  point  are  expressed  in 
terms  of  the  position,  velooity,  and  aooeleration  of 
the  particle  at  a  time  which  is  earlier  than  the 
observation  time,  by  an  amount  sufficient  to  allow. for 
propagation  at  the  velocity  o  from  the  partlole  to  the 

point  of  observation. 

Si  hi  i far 

©a e-same  equations  can  be  used  for  the  corresponding 
advanced  solution,  if  the  partiole’s  position,' velocity, 
and  acceleration  are  taken  from  a  later  point  along  its 
trajectory,  where  the  particle  time,  t^,  and  the  observation 
time,  t,  are  related  by  - 

tx  »  t  +  s/c  ,  *  (1.1) 

ifith-  s  '•  the  propagation  distance.  In  this  case,  the 
radiated  fields  are  received  before  they  have  been  emitted.. 


Pago  iai 


It  is  obvious,  of  course,  that  the  advanced  solution, 
involving  (1*1) »  is  to  be  excluded  on  physical  grounds, 
and  that  only  the  retarded  solution  is  to  be  considered 
as  satisfying  oausality  requirements.  However,  both 
are  exact  solutions  of  Maxwell's  equations,  so  that  the 
selection  of  ene  over  the  other  must  be  made  on  the  basis 

of  an  auxiliary  condition  or  requirement,  which  is  not 

or  implioltly 

explicitly/ contained  in  Maxwell's  equations  themselves. 

Without  this  auxiliary  condition,  a  general  solution  for 
the  electromagnetic  field  associated  with  the  moving 
charged  particle  would  bo  a  superposition  of  the  advanoed 
and  retarded  solutions,  with  an  adjustable  parameter 
specifying  the  proportion:'  in  which  these  two  solutions 
entered  the.  general  solution. 

In  oases  where  exact  solutions  of  Maxwell's  equations 
are  not  available,  numerical  methods  can  be  used,  based 
on  a  finlte-differenoe  analog  of  the  differential  equations. 
However,  the  problem  of  oausality  will  arise  here  again. 

Since  the  differential  equations  do  not  distinguish  between 
advanoed  and  retarded  solutions,  the  flnite-dlfforenoe 
analogs  of  these  equations  will  be  unable  to  make  any 
such  distinction.  As  in' the  oase  of  the  differential  equations, 
an  auxiliary  condition  must  be  imposed,  to  select  the 
retarded  solution. 
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There  is  a  olass  of  dlffioult  but  important  problems, 

in  which  an  imposed  primary  current  of  transient  form 

a 

moves  through  and  near/.regioni:  of  transient  conductivity. 

a  secondary 

The  fields  from  the  primary  our rent  induoe/current;:  in 

>■ 

the  region:  of  finite  conductivity,  and  this-,  secondary 

ourrent  makes'-’  its  own  contribution  to  the  looal 

fields*  v  ••  Henoo  the  secondary  ourrent  -  will  itself L  vs 

of  its  o^m  radiation. 

be  modified  as  a  oonsequence/  An  integral  equation  for 

»» 

the  secondary  ourrent  has  then  to  be  solved.  However, 
all  of  the  field  contributions,  made  by  primary  ourrent 
and  by  secondary  current,  oust  be  retarded  fields,  if  the 

causality  requirement  is  to  be  satisfied. 

in  the  general  case  expressed  as 
While  a  problem  of  this  form  1  s/a  four-dimens lonal.,  ..  . 

space-time  integral  equation,  there  are  certain  special 

symmetry  conditions,  of  Importance  in  practical  examples, 

which  permit  the  problem  to  be  separated  in  polar  coordinates, 

and  reduced  to  a  two-dimensional  problem  involving  only 

a  radial  distance  and  the  time  variable.  Unfortunately, 

further 

the  reduction  to  two  dimensions/obscures  the  causality 
•condition,  if  the  set  of  differential,  equations  is  transformed 
directly  and  then  replaced  by  the  corresponding  set  of 
f inite-dlfferenoe  equations.  An  attempt  to  solve  these 
equations  numerically  leads  to  a  serious  violation  of 
oausality  requirements,  as  will  be  shown  later. 


i 

t 
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A  finite-step  procedure  whioh  utilizes  a  set  of 

lattice  points  in  the  (r,t)-plane,  but  which  satisfies 

oausallty  requirements  at  every  step,  is  also  described 

here.  This  is  a  lattloo  version,  or  finite-difference 

version,  of  an  integral  method  developed  for  the  solution  ' 

«  differential  raethoc 

of  electromagnetic  transient  problems.  This  method,  like  the/ 

makes  use  of  an  expansion  in  veotor  spherical  harmonics,  for 

- 

the  source  ourrents  and  the  field  vectors.. 

The  lattloe  version  of  the  integral  method  is  an 

Iterative  solution  of  the  integral  equation  for  the 

•eoondary  current; \  The  iteration  progresses  outvrard  in 

radial  distance,  and  forward  in  time,  within  the  problem 

domain  in  the  (r,t)-plane.  At  each  step  the  local  electric 

field. is  computed,  as  a  summation  of  the  contributions  of 

primary  and  seoondary  ourrents  previously  computed,  with 

retardation  fully  accounted  for.  The  product  of  the  local 

eleotrlo  field  and  the  local  conductivity  function  (a  knovm 

density 

quantity)  is  then  equal  to  the  seoondary  current/at  that 
lattice  point.  This  seoondary  ourrent  is  then  added  to 
the  primary  ourrent  at  that  lattice  point  (a  known  quantity). 

The  iteration  ’now  progresses  to  the  next  point  in  the  lattice, 

.  and  so  on,  until  the  full  distribution  of  the  seoondary 
ourrent,  over  the  (r,t)-plane,  has  been  determined. 


i 


2.  SECONDARY  CURRENT  C  .  : 

The  oausallty  problem  is  illustrated  most  clearly 
if  the  transient  conductivity  function  does  not  depend 
on  angle,  tut  only  on  radial  distance  and  time* 

*  -  *(r,t)  .  .(C.l) 

It  will  further  be  assumed  that  the  secondary  current  is 
strictly  the  produot  of  the  local  eleotrio  field  and  the 
Icoal  oonduotivity,  independent  of  the  magnitude  of  the 
local  magnetic  field.  In  other  words,  the  Hall  coefficient 
for  the  problem  region  will  be  assumed  to  be  zero. 

Under  these  simplifying  conditions,  the  current 

density  can  be  written  as 
% 

J  =  X  +'  o(r,t)  E  ,  (2.2) 

where  K,  Is  the  primary  current.  For  the  particular 
problem  toward  which  this  analytical  work  has  been  directed 
the  calculation  of  the  electromagnetic  radiation  associated 
with  a  nuclear  detonation  in  the  lower  atmosphere^,  the 

4 

primary  current  JK,  oonslsts  of  Compton  electrons  ejected 
from  air  molecules  by  gamma-rays  from  the  nuolear  reaotion. 
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The  conductivity  cf(r,t)  is  due  to  the  momentary  presence 
in  the  air  of  free  electrons,  released  In  the  Ionizing 
collisions  which  slow  these  Compton  electrons,  plus  the 
later  presence  of  heavier  ions ,  such  as  the  OjT  ions  formed 
when  these  free  electrons  become  attached  to  oxygen  molecules 
in  the  air,  plus  the  positive  ions  left  behind  when  the 
Compton  electrons  were  ejected. 

While  some  bending  of  the  trajectories  cf  these 
oharged  particles  is  to  be  expected,  in  the  transient 
magnetic  field  th  t  is  generated  by  their  motion,  it  will 
be  assumed  (as  noted  above)  that  the  secondary  current 
component  generated  by  this  bending  is  negligible  in 
comparison  with  the  current  J  given  in  (2,2),  jf  Tne 
advantage  of  these  simplifying  assumptions  is  that  there 
results  a  clean  separation  of  the  problem  into  a  sequence 
of  modes  which  are  not  coupled  together.  The  integral 
equation  for  the  secondary  current  separates  into  a 
set  of  uncoupled  Integral  equations,  one  for  each  mode. 
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3.  SEPARATION  INTO  MODES 

When  the  conductivity  function  has  the  form  (2.1), 
implying  that  the  physical  process  to  be  analyzed  has 
a  center  of  at  least  partial  symmetry,  it  is  appropriate 
to  make  an  expansion  in  spherical  coordinates  about  this  ' 

*  L 

central  point.  In  the  example  under  consideration  , 

the  central  point  is  the  location  of  the  nuclear  detonation , 

•through  this 

which  produces  the  primary  current  and/the  transient 
conductivity. 

The  electric  and  magnetic  fields,  E(r,t)  and  H(r,t), 
will-  be  expanded  in  terms  of  the  vector  spherical  harmonics,-* 

and  Sn^9*^*  The  ansular  dependence . is  . 
contained  within  these  harmonics,  while  the  expansion 
coefficients  contain  the  dependence  upon  radial  distance,  r, 
and  time,  t.  The  electric  field,  in  particular,  will  have 
the  expansion 

co  +n 


E(r,t)  -  £  £[ErVr,t)  +.  &<•■*> 

+  E“>n(r,t  (3.1) 


n=0  m=-n 


1 
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The  nfagnetic. field  has  an  expansion  of  the  same  form, 
but  with  the  expansion  coefficients  in  (3.1)  replaced  by 
H£n(r,t),  H®jn(r,t),  and  Hc“n(r,t).  Similarly,  the 
source  current,  J,  and  the  primary  current,  K,  will 
have  expansions  in  terms  of  the  same  vector  spherical 
harmonics,  but  it  will  sometimes  be  convenient  to  write 
the  source  current  as  J(r,  ,t-,),  with  an  expansion 

«M  jl 

involving  suoh  terms  as  J^n(rltt^)  £^ei*^l^»  in  order 
that  the  source  point  (r^,t^)  should  be  clearly  distinguished 
from  the  field  point  (r,t). 

When  there  is  rotational  symmetry  about  the  vertical 
axis,  so  that  the  source  current  does  not  depend  on  the 
azimuth  angle  then  the  fields  will  also  be  Independent 

of  fit  and  the  expansions  such  as  (3.1)  will  be  limited  to 
those  terms  for  which  m  *  0,  Furthermore,  if  the  source 
current  consists  of  radial  and  polar  components  only,  with 
no  current  circulating  about  the  vertical  axis,  then  the 
expansions  will  be  limited  still  further.  With  these 
simplifications,  the  electric  field  (3.1 )  is  reduced  toj 

oo 

E(r,ty  .  £°(e,J5)  +  E°n<r,t)  ,  (3 

n*0 
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while  the  magnetic  field  takes  the  still  simpler  form; 


CO 

H(r,t)  «  ^  •  Hc°>n(r,t)  0°(e,p)  , 


(3.3) 


the  term  with  n»0  feeing  omitted  "because  the  vector 
spherical  harmonic  C®  is  identically  zero,  as  is  3®  , 
though  Pq  does  not  vanish. 

The  explicit  forms  for  the  vector  spherical  harmonics 
in  (3.2)  and  (3.3)  are 


*£<•»*> 

Eg<e.P> 

c°(e.p) 


ap  Pn(cos  6)  , 
-4==  Pl(oos  ©)  , 


Vn(n+1)  n 


'n(n+l)  n 


Pi  (cos  e)  , 


(3.4) 


(3.5) 


(3.6) 


where  Pn(cos  6)  are  the  Legendre  polynomials,  and 


pl(oos©)  =  -  —  Jn(00S  S ) 


sln©v7 - r-r-  P„( cos©)  , 

d(cos  e)  n 


(3.7) 


The  unit  vectors  ar,  aQ,  a^  oan  fee  expressed  in  terms 
•of  the  rectangular-coordinate  unit  vectors,  i,  J,  k,  through; 
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5© 

a  01 

-P 


1  sin  6  cos  p 
i  cos©  cos p 
-  i  sin  p  + 


+  J  sin©  slnp  +  k  oos©  , 
+  ^  cos  ©  sinp  -  k  sin©  , 
1  oos  p  . 

«M  9 


(3.8) 


The  Legendre  polynomials  satisfy  the  second-order 
differential  equation 


d* 

d©2 


Pn(cos  ©) 


+  cote 


d© 


Pn(oos  ©)  +  n(n+l)  Pn(cos  ©) 


It  will  be  assumed  that  the  medium  is  the  vacuum,  and 
that  all  dielectric  and  conduction  effects  are  interpreted 
through  secondary  currents  that  will  be  explicitly  computed 
during  the  calculations.  Furthermore,  it  will  be  assumed 
that  the  dielectric  effects  can  be  neglected  in  this  analysis , 
and  that  the  conduction  effects  can  be  incorporated  into  the 
source-current  distribution  through  Eq,  (2*2). 


With  these  assumptions,  Maxwell's  equations  take  the  form 
S*S  =  %  +  5!  +  60^E  (3.10a) 

•  (3.10b) 


in  KKS  units.  When  the  primary  current  K  is  given  an 
expansion  of  the  form  (3.2),  and  when  3  and  H  are  represented 
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by  (3.2)  and  (3.3) *  substitution  into  (3.10)  leads  to 

2 

the  three  scalar  equations  ; 

•.|&[*  HC?n<r*fc>]  "  +  *»*<*•*' 

+  £o  &>*.„<’■*’  ’  (3*11) 

iyiltn+l)  Hc°in(r,t)  =  fc£n<r,t>  +  Er,n(r,t) 

+  €.  ^  Er,n<^,t,  ■  0.12) 

(r,t) 
(3 


O  Qt 


HC°,n<r’t’  “  F 


[r  EB°n(r>t,J  -  r^n(n+1)  Sr,n 


When  successive  integer  values  of  n  are  substituted,  feach 
mode  is  then  described  by  a  set  of  three  coupled  first-order 
differential  equations  relating  the  unknown  field  components, 
Er°.n  •  EB°,n  •  HC.n  •  to  the  (presumed  to  be  known)  primary 
current  components ,  ’  K^?  n  and  Kg  n  ,  and  the  (presumed  to 
be  known)  conductivity  function,  cr(r,t). 

Equations  (3.11-13)  are  a  transformed  version  of  Maxwell's 
equations  (3.10),  and  therefore  admit  both  retarded  and 
advanced  solutions,  or  arbitrary  linear  combinations  of 
retarded  and  advanced  solutions.  The  physical  causality 
requirement .is  not  contained  in  (3.11-13) $  and  must  be 

specified  through  an  auxiliary  condition,  if  mathematical 
solutions  of  (3.11-13)  are  to  represent  actual  phenomena. 
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4.  RETARDED  HERTZ  VECTOR 

A  solution  to  Maxwell’s  equations,  and  in  particular 
a  solution  to  the  soalar  equations  (3.11-13),  oan  be  based 
on  the  retarded  Hertz  vector.  In  this  way  the  limitation 
to  the  retarded  solution,  with  the  complete  exclusion  of  • 
the  advanced  solution,  i^  fully  enforced. 

In  terms  of  the  Hertz  vector,  TT  »  the  electric  and 

A) 

magnetic  field  vectors  are  expressed  through^: 

5  -  2Y‘TT  -  4  -^TT  .  <*.i) 

5  '  •  <*-2> 

The  retarded  Hertz  vector  is  a  function  of  the  actual 
source-current  distribution,  J,  as  given  in  (2,2)  for 
the  special  symmetry  which  is  being  considered  here,^ 

In  the  general  case,  with  the  source- current  density 
J(ri,ti)  given  the  full  expansion,  involving  the  vector 
spherical  harmonios  P®,  bP,  and  C®,  the  retarded  Hertz 

a*n  mu  Mil 

vector'  separates  into  modes  corresponding  to  the  modes 
for  the  source  current.  For  each  mode,  the  Hertz-vector 
expansion  coefficients  are  given  in  terms  of  the 
source- current  expansion  coefficients  by  the  matrix  equation 
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r,=°o 


HV“n<r,t)  =  ?!;  J  (VrlJrl^l<r'rl-rwl1  Jj£,n<rl’h>dtl 


rl=0 


tl=to 


rx=r  tj=t-(r-ri)/o 

JcVr>driO<>.ri.r-ri)  J  ^(r^t^dt^ 


rl=° 


r-=oo 


ti=t0 


-  tj=t- (r^-r )/ o 

J(Vr)drlGa^(r.r1,r1-r)  J  J^Ov^Jdtj 


rl=r 


rlP 


h=1o 


26 


-I/--  t^t- (r-rx  )/c 

j  (rj/r  )dr^  J  ^(r.r^et-otj.)  dtx 

t^t-fr+r^J/c 


rl=0 


ri?"  hr1'  !ri'r  )/c 


‘lr  1*  "  v 

Jtr/rjdrj  J  O^tr.rj.et-o^)  dtj 

tj=t-(r+r^)/c 


rl=r 


(4.3) 


The  Greek  indices,  <J  and  X,  run  through  the  three  values: 
r,  A  repeated  Greek  index  indicates  summation  over  these 

three  .values.  Thus  (4.3)  represents  three  equations,  according 
as  represents  r,  B,  or  C. 
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In  the  Integrals  in  (4.3 )»  the  time  tG  is  chosen  to 
precede  the  Initiation  of  the  transient  current  flow,  whioh 
is  responsible  for  the  generation  of  the  field  distribution 
which  is  incorporated  into  the  Hertz  vector  (4.3)*  Any 
pre-existing  static  electric  field  can  be  expressed  as 
an  earlier  flow  of  current,  producing  the  distribution  of 
electric  charge  which  is  to  be  associated  with  this  initial 
electrostatic  charge .distribution.  Thus  if  the  initial 
moment,  tQ,  is  pushed  back  far  enough  in  the  time  scale, 
the  Hertz  vector  (4.3)  will  incorporate  all  transient 
effects  and  any  pre-existing  static  effects. 

*  to- 

Once  the  retarded  Hertz  vector  (4.3)  has  been  expressed 

in  terms  of  a  particular  source- current  distribution, 

the  corresponding  electric  and  magnetic  field  vectors  can- 

be  calculated  directly  through  (4.1-2),  and  it  is  found 

that  the  mode  separation  is  retained.  For  the  particular 

symmetry  described  by  (3«11-13)»  the  three  nonvanishing 

electromagnetic  field  components  are  all  found  to  have 

explicit  representations  of  the  same  general  form  as  (^.3). 

Eaoh  of  the  components  EjP n(r,t) ,  EB,n^r*^*  HQ>n^r*^  * 

is  expressed  as  a  sum  of  several  integrals  containing 
e  0  with  the  integration 

Jr  n^l’H^  or  JB  n^rl*^l^  intesrand, /taken  over 

*  3 

the  causally  accessible  region  of  the  (r^ jt^) -plane. 
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5.  RETARDED  FIELDS  FOR  FIRST  MODE 

The  mode  for  which  n*0,  which  will  be  called 

the 'zeroth-  •v.r  mode,'  Involves  only  the  radial  component 
field 

of  the  eleotrlc/and  the  radial  component  of  the  ourrent 
density.  The/ equation  that  relates  these  radial  components  is 

*,0,0(r.t)  -  -  ^  /  Jr,o<r*V  dtl  *  <5.U 

°  V*o 

whioh  Is  equivalent  to  Gauss’s  law/  For  this  mode  the 

present 

effects  of  retardation  are  /.;•  In  the  regarded  Hercz  vector, 
but  drop  out- when  the  eleotrio  field  (5.1)  Is  computed 

from  the  Hertz  vector, 

•  • 

However,  retardation  effects  remain  in  the  field 

components  that  are  computed  for  the  mode  with  n*l, 
whioh  will  be  oalled  the  ’first  mode,*  The  distinction 
between  oausal  and  non-oausal  solutions  of  Maxwell’s 
equations  can  thus  be  Illustrated  by  the  fields 
associated  with  this  first  mode.  It  will  be  assumed, 
as  stated  earlier,  that  the  oirculatory  current  component* 

^<?,1  ’  *s  zero»  80  only  the  field  components  Ej?  1  »  EB  1  * 

and  are  generated.  The  retarded,  fields,  obtained  from 

theretarded  Hertz  vector,  will 'then  be  written  explicitly, 
with  the  aid  of  the  abbreviation; 


Page  136 


r, -r  t,«t-(r-r,  )/o  r,»co  -r)/o 


(If) 


XL 


l?  'l 

J  J 

“  v  »  »  a  a 


(5.2a) 


rl*0  t1*t-(r+r1)/c  r^«r  t^t-fr+r^)/© 


whioh  oan  also  be  written  In  the  equivalent  form  (but  with 
the  order  of  integration  reversed); 


ti*t-r/o  r^r+oft-t^)  tj»t  r^r+oft-t^) 

(//)  -  j  T  +  j  T  .  (5. 

tjB-oo  r^= -r+oft-t^)  t^t-r/o  r^r-cCt-t^) 

In  (5.2b)  the  limit  t^ss-co  will  ordinarily  be  replaced  by 
tj-t©.  Where*  t©  is  a  time  preceding  any  of  the  transient 
eurront  flow. 

The  radial  component  of  (4*1),  for  this  first  mode,  gives 
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v! 
3  €c 


ri?'°  ti^t-(rx-r5/c 


J  J  dtidri 

rx«r  tx-t0 


“  tV  ^  [(ri2/r3)  +  (1/rlO  Jr,l(rl*V  dtidrl 

-  (//)  [2(rx?/r3)  -  (l/rx)]  JB,l(rl»V  dtl  dri 

♦  J^T  (//)  (1/ r3rx )  ^3(r2+r12  )c(t-t1>  -  c3  (t-t^l  j^Cr^t-^  dtl  dr 

-  g-g-  (//■)  (l/r3^)  |^3(r2-r12  -  o3(t-t1)3JJ®>:L(r1,t1)  4^  dr 

(5.3; 


The  component  of  Bq,  (4,1)  which  is  in  the  polar  direction 

**  V 

gives  the  result,  for  this  first  mode; 

EB,l(r,t>  =  Fl~  J  (rx/or)' JB>l{r1»t-[p*riyo)  dPj 

°  rl=°. 


(r^/cr)  J 


0 

B,1 


(r^t-fr-rj/c)  drx 


(rj/cr)  jb  1(r1»t-[r1-r]/c)  dr. 


1 

2 
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-ij-  j  J  '(*//»’>  [j^.V  +  ^  ^.x'W]  dtl  drX 


t,=t-(r-r,  )/c 


*1=°  VtQ 


r17>m  tl?t-(rX"r,/o 


3  6. 


(l/r 


1>  [v? 


-  <1^1’^’]  dti  dri 


v? 


■  6  6, 


(//)  [(r^/r3)  -  2  (l/r,)]  j£i<*i.V  dtl  drl 


dr. 


+  3T"  (//)  [(rl2/r3)  +  (l/rl>]  JB,l(rl’tl)  dtl  dri 

+  (//)  (l/r3r1'  jjHr^-r^  )c(t-t1'  +  o3(t-t1)3jj£fl<r1»t1)dt1 

~  (//I  (X/r3rx)  r^C^+rj2  )c(t-tx)  +  o3(t-t1)3Jjgjl(r1,t1)dt1dr1 

(5» 


12  6, 


The  nonvanishing  component  of  Eq.  (4*2) »  which  is  in  the 
azimuthal  direction’ a ^  ,  gives  the  corresponding  equation 


I  \ 


for  the  nonvanishing  component  of  the  magnetic  field; 


i 
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rlssC0 

HC,l(r,t)  *  “  |  J  (rl/r)  JB,l(rl’Mr+rll/o)  ^1 

r^O 

rjrr 

+  |  /  (rj/r)  [r_ri]/cl  drx 

0 

Ty=CO 

-  |  J  (rx/r)  •JB,i<r1,t-[r1-r]/o)  drx 
+  (//)  (l/r2^)  ^(i^+r.2  )  -  c^Ct-t^J  dtl  drX 

-  *5“  (//)  (l/r2^)  J^r2-^2)  -  o2(t~t1)2J  JQjl(ri>ti)  drx  . 

(5.5) 

It  can  be  verified  through  explicit  substitution  that  the 
three  field  components  for  this  first  mode,  as  given  in 
(5.3-5) »  satisfy  the  reduced  equations  (3.11-13)*  if  the 
source  current  is  given  by  (2,2), 


( 


/ 
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6.  NUMERICAL  SOLUTION  OP  INTEGRAL  EQUATION’ 

•  -9 

The  integral  equation  tfo  be  solved  is  actually  a 
coupled  pair  of  equations  for  the  current-density 

=  +■  tt (r, t)  EjP^fr.ty  ,  (6.1a) 

=  Katl(r,t)  .+  ,<»(r,t)  ,  (6.1b) 

together  with  the  equations  (5.3)  and  (5.4)  which  give 
the  electric  field  components  in  terms  of  integrations 
over  the  current-density .components, 

4 

In  the  problem  of  direct  Interest  ,  the  primary  current  K 
and  the  transient  conductivity  function  are  both  initiated 
by  the  nuclear  reaction,  and  are  therefore  zero  before  the 
time  of  the  detonation.  At  a  distance  r  from  the  detonation 
point,  K  and  will  be  zero  up  until  a  time  which  is  later 
than  the  detonation  time  by  the  amount  r/c,  the  propagation 
time  for  gamma  rays  moving  radially  outward  from  the  detonation. 
If  the  detonation  occurs  at  t  =  0,  then  the  problem  domain 
in  the  (r,t)-plane  is  limited  to  the  region  to  the  right  of 
a  diagonal  line  through  the  origin,  as  shown  in  Figure  1, 


components, 

J°>x(r,t) 

jfA(r.t) 
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For  the  numerical  solution,  the  problem  domain  is 
subdivided  by  a  diagonal  lattice,  bounded  on  the  left 
by  the  line 

r  =  ct  ,  (6.2a) 

and  bounded  below  by  the  line 

r  =  0 .  (6.2b) 

Field  and  current  component  values  are  calculated  for  the 
points  where  lattice  lines  cross.  Integrations  over 
previously- calculated  current  values  are  expressed  as 
summations  over  the  small  square  areas  in  the  lattice. 

The  lattice-line  spacing  is  made  sufficiently  fine 

so  that  the  error  introduced  through  the  discreteness 

is  small.  In  practice,  the  line  spacings  can  differ  in 

different  portions  of  the  problem  domain,  and.  the  small 

areas  can  be  rectangular  Instead  of  square,  If  necessary 

to  make  efficient  use  of  computing-machine  time  and  storage. 

At  the  domain  boundary  (6.2a),  in  particular,  the  lattice 

spacing  should  be  small  enough  so  that,  in  the  problem 
k 

being  analyzed  ,  the  secondary  current,  d  3,  is  still  very 
small  in  comparison  with  the  primary  current,  K.  The 
initial  calculations  of  the  field  components  can  then  be 
based  on  the  known  primary  current.  As  the  computation 
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progresses  from  one  lattice  point  to  the  next,  and  moves 
to  regions  that  are  not  close  to  the  boundary  (6.2a), 

4 

the  secondary  current  will  play  an  Increasingly  Important 

role. 

The  numerical  solution  also  requires  Initial 
ourrent-denslty  values  along  the  inner  boundary  (6* 2b), 

h, 

In  the  problem  being -considered  ,  this  boundary  refers 
to  the  immediate  neighborhood  of  the  nuclear  detonation, 
for  times  subsequent  to  the  detonation  time.  The  asymmetry 
which  leads  to  the  first-mode  portion  of  the  primary  current 
is  the  asymmetry  produced  by  the  presence  of  the  ground 

9 

or  by  the  presence  of  an  atmospheric  air-density  gradient. 

In  practice  this  asymmetry  will  not  enter  until  at  least 
a  short  distance  from  the  detonation  center.  Thus  the 
•first-mode  portion  of  the  primary  current  can  be  set  equal 
to  zero  along  the  inner  boundary  (6.2b),  in  this  application 
of  the  theory.  For  other  applications,  physical  considerations 
or  a  simple  analytical  model  will  generally  give  an  adequate 

basis  for  setting  the  initial  current-density  values  along 

\ 

the  boundary  (6.2b). 

CflTv  'pr'OC enfcl 

The  numerical  iteration  proceeds  diagonally  upward 
to  the  right,  in  f?g.  1,  alons  lattice  lines  parallel  to 
the  diagonal  boundary  (6.2a),  One  step  in  the  iteration 
is  indicated  in  Figure  2.  Here  the  current  components  (6.1) 
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are  being  computed  at  the  lattice  point  P,  with  the 

Integrations  (5*3)  and  (5*4)  replaced  by  summations 

over  the  lattice  squares  shown  in  the  figure.  Each 

square  Is  treated  as  though  It  were  concentrated  into 

a  point  souroe  having  the  location  of  its  central 

point,  shorn  here  as  a  dot  in  the  center  of  the  square. 

The  current  associated  with  the  square  is  the  average 

of  the  current  values  at  its  oorners,  previously  computed, 

(The  triangles  at  the  bottom  can  be  included  separately. ) 

In  the  case  of  the  square  which  has  the  point  P  at 

its  right-hancl.  corner,  the  unlcnown  current  at  P  oan  be 
approximated  through  extrapolation  from  the  three  currents 
at  the  other  oorners.  Alternatively,  the  current  at  P 
oan  be  included  formally  in  the  summation  (to  which  the 
integration  has  been  reduced  by.  this  approximate  method) , 
then  transferred  algebraically  to  the  left-hand  side  of  (6.1), 
and  solved  for  as  a  part  of  the  numerical  solution  for  the 
current -density  components  at  P. 

Either  alternative  will  give  the  current  at  P,  and 
the  iteration  can  then  step  to  the  next  point  on  the 
lattice.  The  final  result  will  be  a  set  of  current-component 
values  for  all  the  lattice  points.  If  desired,  the  program 
oan  then  be  re-run,  with  these  values  substituted  when  the 
current  at  P  is  needed  for  the  integral  over  the  square 
eontainlns  the  point  P  as  its  right-hand  corner.  If  the 
lattice  is  fine  enough,  the  change  in  current  values,  resulting 
-from  thi?  re-run,  should  be  small. 
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The  numerical  solution  which  has  been  described 
will  give  the  current  values  at  the  lattice  points, 
when  the  primary  current  K  and  the  conductivity  c f  are 
inserted,  along  with  the  boundary  values  for  the 
ourrent  on  the  inner  boundary  (6.2b).  The  transient 
eleotrio  and  magnetic  field  components  will  also  be 
given,  through  integration  of  this  current  distribution, 
using  (5.3-5).  The  result  will  be  consistent  with 
causality  requirements,  provided  the  lattice  is 
sufficiently  fine.  In  practice,  if  the  solution  is 
not  changed  significantly  when  the  lattice  is  made 
still  finer,  then  it  was  already  sufficiently  fine. 

While  this-  illustration  has  referred  to  the  first 
mode,  the  mode  with  n«l,  the  same  general  considerations 
apply  to  the  solution  of  the  integral  equation  for  any 
higher  mode,  as  long  as  the  modes  are  separable.  If 
the  conductivity  function  is  not  given  by  d(r,t),  but 
by  a  function  which  depends  upon  ©  as  well,  then  the 
modes  will  be  coupled  together  and  (6.1)  will  be  replaced 
by  an  extended  matrix  equation,  coupling  together  the' 
current  compcnents  for  the  different  modes.  In  this  case, 
a  generalized  version  of  the  iteration  described  here  can 
still  be  used,  without  violating  the  causality  requirements. 
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7.  NON- CAUSAL  SOLUTION 

2  1 

For  comparison,  a  non-causal  solution  of  Maxwell's 
equations,  in  the  reduced  form  given  in  (3.11-13)*  will 
be  desoribed  briefly.  In  this  non-causal  solution, 
the  sane  pr obi era  domain  of  Fig.  1  has  been  used,  but 
the  lattioe  has  the  orientation  shown  here  in  Figure  3. 

The  iteration  progresses  vertically  upward  from  the 

% 

boundary  (6.2b),  and  requires  initial  values  along 

this  boundary. and  also  along  and  near  the  boundary  (6.2a). 

An  iteration  step  is  shown  in  Figure  4.  The  field 
components  are  caloulated  at  the  point  P,  based  on  the 
primary  current  and  the  conductivity  at  the  point  X 
and  the  previously-calculated  field  components  at  the 
three  points  Q.  A  finite-difference  version  of  the 
equations  (3.11-13)  is  used  in  this  calculation.  The 
result  Is  a  numerical  solution  of  Maxwell's  equations, 
for  each  step,  but  each  step  gives  a  mixing  of  advanced 
and  retarded  solutions  since  the  causality  requirements 
•have  not  been  Imposed  in  any  way.  Thus  the  result  is 
a  mathematical  solution  of  Maxwell ' s  equations  vrhich  has 
a  doubtful  relationship  to  physical  phenomena, 

f 
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It  can  be  seen  from  Fig.  4  that  the  introduction  of 
an-  intensification  of  the  primary  current  at  the  point  X 
that  is  shown  will  have  effects  on  point  P,  as  shown ,  and.  hence 
also  on  other  points  P  which  lie  directly  above  this  one 
in  the  figure.  Thus  a  current  change  at  an  inner  radius 
will  lead  to  field  changes  at  an  outer  radius  which  are  -, 
essentially  simultaneous,  and  must  have  traveled  faster  than  c, 
These-effects  are  physically  inadmissible  ,  yet  they  do 
satisfy  the  finite-difference  form  of  Maxwell's  equations 
which  has  been  used  in  the  iteration  step  of  Fig.  4. 

The  solution,  of  course,  is  not  to  change  to  another 
form  of  iteration  step  which  just  translates  Eqs-.  (3.11-13) 
from  differential  form  to  difference  form  in  some  other 
way.  The  solution,  obviously,  .to  the  problem  of  solving 
Maxwell's  equations  numerically  in  the  (r,t)-plane, 
is  to  impose  the  causality  requirements  at  an  earlier 
stage,  when  the  equations  are  in  four-dimensional  space-time, 
and  to  select  the  retarded  solution  at  the  beginning* 

The  reduction  to  the  (r,t)-plane  then  leads  to  the 
retarded  Hertz  vector  (4,3)  and  the  retarded  fields 
illustrated  in  (5.3-5). 
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8.  CONCLUSIONS 

In  the  general  case,  a* numerical  solution  of 
•Maxwell's  equations  will  give  a  superposition  of 
electromagnetic  phenomena  which  can  occur  in  a 
physical  situation, and  other  phenomena  which  cannot 
occur  physically.  If  the  mathematical  solution  is 
to  describe  a  physical  process,  then  the  analysis 
must  be  set  up  in  such  a  way  that  only  the  retarded-field 
solution  of  Maxwell's  equations  is  admitted,  while  the 
advanced  fields  are  excluded  from  the  beginning. 

In  problems  where  there  is  a  center  of  symmetry  or 
of  partial  symmetry,  a  description' of  the  problem  in  terms 
of  spherical  polar  coordinates  is  appropriate.  In  certain 
problems  an  expansion  in  vector  spherical  harmonics 
provides  a  separation  into  modes,  and  a  reduction  of  the 
four-dimensional  problem  into  a  set  of  t*jro- dimensional 
problems,  one  for  each  mode.  In  the  reduced  problem,  the 
variables  are  the  radial  distance  r  and  the  time  t. 

If  retardation  is  ignored,  then  a  non-causal  solution 
In  the  (r,t) -plane  can  be  obtained  through  a  point-to-point 
numerical  integration  of  a  finite-difference  version  of 
Maxwell's  equations,  as  they  appear  when  reduced  to  two 
dimensions  by  the  separation  Into  vector-spherical -harmonic 
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modes.  This  solution  will  in  the  general  case  violate 
the  causality  requirement  and  must  therefore  be  ruled 
out  on  physical  grounds. 

When  retardation  is  incorporated  at  the  beginning 
of  the  analysis,  before  the  separation  of  variables, 
then  the  physical  causality  requirement  will  be  satisfied.' 
The  reduction  to  the  (r,t) -plane  can  still  be  carried  out, 
but  in  this  case  the  result  is  an  integral  equation  to 
be  solved  numerically,  rather  than  a  differential  equation 
The  solution  of  the  integral  equation  will  provide  one 
of  the  possible  solutions  of  the  differential  equation, 
but  will  exclude  the  non-physical  solutions. 

The-  integral  equation  can  be  solved  numerically 
with  the  aid  of  a  discrete  lattice  in  the  (r, t)-plane. 

Each  iteration  step  then  involves  a  summation  over  all 
the  causally-accesslble  lattice  positions  in  the  region 
covered  by  previous  iteration  steps.  The  iteration 
procedure  resembles . in  a  mechanical  way  the  iteration 
used  in  the  non-causal  finite-difference  solution,  so  that 
the  conversion  of  a  computing-machine  program  from  the 
,  non-causal  solution  to  the  causal  solution  can  readily 
be  carried  out. 
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